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Abstract 



Gross-Neveu type models with a finite number of fermion flavours are studied 
on a two-dimensional Euclidean space-time lattice. The models are asymp- 
totically free and are invariant under a chiral symmetry. These similarities to 
QCD make them perfect benchmark systems for fermion actions used in large 
scale lattice QCD computations. The Schrodinger functional for the Gross- 
Neveu models is deflned for both, Wilson and Ginsparg- Wilson fermions, and 
shown to be renormalisable in 1-loop lattice perturbation theory. 

In two dimensions four fermion interactions of the Gross-Neveu models 
have dimensionless coupling constants. The symmetry properties of the four 
fermion interaction terms and the relations among them are discussed. For 
Wilson fermions chiral symmetry is explicitly broken and additional terms 
must be included in the action. Chiral symmetry is restored up to cut-off 
effects by tuning the bare mass and one of the couplings. The critical mass 
and the symmetry restoring coupling are computed to second order in lattice 
perturbation theory. 

This result is used in the 1-loop computation of the renormalised cou- 
plings and the associated beta-functions. The renormalised couplings are 
defined in terms of suitable boundary-to-boundary correlation functions. In 
the computation the known first order coefficients of the beta-functions are 
reproduced. One of the couplings is found to have a vanishing beta-function. 
The calculation is repeated for the recently proposed Schrodinger functional 
with exact chiral symmetry, i.e. Ginsparg- Wilson fermions. The renormal- 
isation pattern is found to be the same as in the Wilson case. Using the 
regularisation dependent finite part of the renormalised couplings, the ratio 
of the Lambda-parameters is computed. 
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Zusammenfassung 



In dieser Arbeit werden Gross-Neveu Modelle mit einer endlichen Anzahl 
von Fermiontypen auf einem zweidimensionalen Euklidischen Raumzeitgit- 
ter betrachtet. Modelle dieses Typs sind asymptotisch frei und invariant un- 
ter einer chiralen Symmetrie. Aufgrund dieser Gemeinsamkeiten mit QCD 
sind sie sehr giit geeignet als Testumgebungen fur Fermionwirkungen die 
in grofsangelegten Gitter-QCD-Rechnungen benutzt werden. Das Schrodin- 
ger Funktional fiir die Gross-Neveu Modelle wird definiert fiir Wilson und 
Ginsparg- Wilson Fermionen. In l-Schleifenstorungstheorie wird seine Renor- 
mierbarkeit gezeigt. 

Die Vier-Fermionwechselwirkungen der Gross-Neveu Modelle habe di- 
mensionslose Kopplungskonstanten in zwei Dimensionen. Die Symmetrieei- 
genschaften der Vier-Fermionwechselwirkungen und deren Beziehungen un- 
tereinander werden diskutiert. Im Fall von Wilson Fermionen ist die chirale 
Symmetrie explizit gebrochen und zusatzliche Terme miissen in die Wirkung 
aufgenommen werden. Die chirale Symmetrie wird durch das Einstellen der 
nackten Masse und einer der Kopplungen bis auf Cut-off-Elfekte wiederher- 
gestellt. Die kritische Masse und die symmetriewiederherstellende Kopplung 
werden bis zur zweiten Ordnung in Gitterstorungstheorie berechnet. 

Dieses Resultat wird in der 1-Schleifenberechnung der renormierten Kopp- 
lungen und der zugehorigen Betafunktionen benutzt. Die renormierten Kopp- 
lungen werden definiert mit Hilfe von geeignete Rand-Rand-Korrelatoren. 
Die Rechnung reproduziert die bekannten fiihrenden Koeffizienten der Beta- 
funktionen. Eine der Kopplungen hat eine verschwindende Betafunktion. Die 
Rechnung wird mit dem vor kurzem vorgeschlagenen Schrodinger Funktional 
mit exakter chiraler Symmetrie, also Ginsparg Wilson Fermionen, wiederholt. 
Es werden die gleichen Divergenzen gefimden, wie im Fall von Wilson Fermio- 
nen. Unter Benutzung des regularisierungsabhangigen, endlichen Teils der 
renormierten Kopplungen werden die Verhaltnisse der Lambda-Parameter 
bestimmt. 



Schlagworter: 

Chirales Gross-Neveu Modell, Ginsparg- Wilson Fermionen, Schrodinger 
Funktional, Gitterstorungstheorie 
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Chapter 1 
Introduction 



The theory of strong interactions, quantum chromodynamics (QCD), is a 
gauge theory with few parameters and is though assumed to describe many 
phenomenons, such as the mass spectrum of the hadrons and scattering pro- 
cesses involving quarks and gluons. Experimental data is available from the 
low energy regime of the light meson masses to the high energy regime of 
hadron-hadron scattering 



At high energies the fundamental degrees of freedom, the quarks and the 
gluons, are only weakly coupled. In such a situation a perturbative treat- 
ment, where the interactions of the quarks and gluons are small corrections, 
is justified. Indeed, perturbative calculations successfully describe the data, 
for example, of deep inelastic hadron-lepton scattering. At low energies the 
coupling becomes large and the interactions are not small but dominant. 
Therefore the perturbative approach is not applicable in this regime. In ad- 
dition, the relevant degrees of freedom are no longer the fundamental quarks 
and gluons, but the lightest bound states (the light mesons, i.e. pions). Chi- 
ral perturbation theory has been developed to accommodate this. But as 
an effective theory it has to be matched to the experiments, thus losing the 
appeal of a first principle computation. 

The lattice discretisation of quantum field theories is a powerful method 
that was first applied to QCD long ago by Wilson Q. Since then lattice QCD 
has proved to allow for non-perturbative calculations from first principle and 
to connect the low and high energy regime (see Q for an example). First of 
all the lattice regulates the theory, in which the inverse lattice spacing 1/a 
serves as a sharp momentum cut-off and thus renders the theory ultravio- 
let finite. In the infrared a non-vanishing mass or a finite volume scheme 
with specific boundary conditions like the Schrodinger functional (SF) pro- 
vides a lower bound on the modes in the theory. (The merits of the SF 
are described in Chapter [H) In this way a quantum field theory on the lat- 
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tice is mathematically well defined without reference to perturbation theory. 
Nevertheless perturbation theory can be used at this point. Being in general 
more complicated than similar studies in the continuum, lattice perturbation 
theory is needed, for example, to translate lattice results into the language of 
continuum renormalisation schemes like MS (minimal subtraction), that are 
mostly used by experimentalists. Furthermore, in the weak coupling regime, 
it serves as guidance and cross check for non-perturbative methods. 

If the metric of the lattice theory is the Euclidean one, the path integral 
representation of quantum field theory is accessible to numerical evaluation 
via Monte Carlo simulations. Today this is the prominent approach to extract 
physics from lattice QCD. Since the days of Wilsons proposal there has been 
substantial progress in the understanding of lattice QCD as well as in the 
algorithms that are used to perform the simulations. Still, the by definition 
limited computational resources are the main obstacle to accomplish the goal 
of producing predictions without any compromise. 

Along the way so called "toy models" were studied. The term refers to 
quantum field theories that are simpler but in some respects similar to QCD. 
They were used, for example, to test new methods [J| or to conjecture the 
phase structure of lattice QCD One reason to consider these simpler 
models is that there might be analytical tools at hand that allow one to solve 
the model exactly. Often another advantage is that Monte Carlo simulations 
of the toy model are much more cheaper. The knowledge and results gained 
in the simpler models can then be used to argue in the more complicated 
theory. Or assumptions and approximations used in lattice QCD without a 
chance to prove their validity there, can be applied and confronted with the 
exact result and/or high precision numerical data in the simpler theory. 

For example, in order to make predictions about the real world, the lat- 
tice discretisation, as any regulator, has to be removed in the end. Numerical 
simulations of lattice QCD are only possible at finite lattice spacing a. In 
practice one computes the observable of interest for a number of lattice spac- 
ings and extrapolates to the continuum limit. If the measured points show a 
significant dependence on the lattice spacing one has to assume a functional 
form to perform the extrapolation. The only known prescription of the lat- 
tice artefacts goes back to the work of Symanzik j^, 0, @, 0|. His conclusions 
for the functional form of lattice artefacts in lattice field theories are based 
on an effective theory and perturbation theory. Using this form in the ex- 
trapolation step of lattice QCD computations introduces a possible source of 
systematic errors in a presumably first principle computation. In the two di- 
mensional and asymptotically free non-linear sigma models these aspects can 
be studied with high precision Monte Carlo simulations and analytic tools 
[lo| . The message for lattice QCD is clear: try to avoid the extrapolation 
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step. This can be achieved partly by eliminating the leading lattice artefacts 
by implementing Symanzik's improvement program 0, @|. 

At the time of writing this thesis a number of collaborations of lattice 
physicists are simulating full lattice QCD with two or three light quarks and 



have presented first results [ill, [iJ, ll3|, llJ, ll5|, [la, ll7| . The main difference be 



tween the approaches followed by these collaborations is in the used fermion 
action. Beside the already mentioned Wilson fermions, there are Wilson 
twisted mass, overlap, staggered fermions and the fixed point action. There 
are also differences in the gauge action and the various kinds of improvement 



18fl 



applied, but here we concentrate on the fermionic part of the action (see 
for an overview). 

Some of these fermion actions are more theoretically sound than others. 
The hope is, based on universality arguments, that in the continuum limit, 
where the correlation length diverges, differences on the scale of the lattice 
spacing are unimportant. Clearly, a numerical test of this presumed agree- 
ment in the continuum limit is desirable. Due to the restricted numerical 
resources of today this is impossible in lattice QCD in the near future. In 
such a situation a two dimensional non-trivial fermionic quantum field theory 
could serve as a benchmark system. If the actions agree in the toy model, 
it would not be a proof for QCD, but some confidence would be gained. If 
they do not agree, it would be clear that there are serious problems and that 
most probably also lattice QCD simulations are affected. 

In this work we study models of self-coupled fermions in two dimensional 
space-time. There are several theories of this kind in two dimensions referred 
to as Gross-Neveu [l9| and Thirring models [20|. We consider here the first 
type. Among them the most similar to QCD is the chiral Gross-Neveu model 
(CGN) with N types (in the following referred to as fiavours) of fermions. 
They are coupled through quartic interaction terms. Since in two dimensions 
the fermion fields have mass dimension 1/2, the corresponding couplings are 
dimensionless. The CGN shares with QCD the features asymptotic freedom 
and a continuous chiral symmetry (in the massless theory). The continuum 



CGN has been studied in perturbation theory up to three loops [2l|, [22 



23l |. Beyond perturbation theory the 5'-matrix and the particle spectrum are 



known to some extent [2J, l25|, l26|, l27|, l28|. Many properties of the model can 



be studied in the limit of infinite many fiavours (large- iV limit) , which is also 
sensitive to the non-perturbative nature of the model. In this limit the model 
is asymptotically free and a fermion mass is dynamically generated [l9|, [29 



On the lattice the model has been studied so far almost exclusively in 



the large- limit [30, iMfl- In this thesis we define the CGN with a finite 
number of fermion fiavours on the lattice. For the fermion action we use 
Wilson's version f3| since it is the most rigorous and theoretically sound 
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one. After the theory has been established in this way, it can be used to 
check other actions. As a first application we analyse a recently proposed 
Dirac operator |32|| that is expected to be compatible with the Schrodinger 
functional boundary conditions and, at the same time, is a solution to the 
Ginsparg- Wilson relation |33| in the bulk of the lattice (up to exponentially 
decreasing tails). A Dirac operator that satisfies this relation has better 
chiral properties than standard Wilson fermions and is thus better suited in 
cases were chiral symmetry plays an important role. Since Ginsparg- Wilson 
fermions are very expensive in terms of computational costs, precise studies in 
two dimensions are very welcome. We define observables suitable for Monte 
Carlo simulations |3j| and compute them in first order lattice perturbation 
theory. This is the first computation with Ginsparg- Wilson fermions in the 
Schrodinger functional beyond the free theory. 

This thesis is organised as follows. In Chapter [2] we give a short review 
of the perturbative renormalisation and discretisation of quantum field the- 
ories. In Chapter [3] aspects of chiral symmetry on the lattice are addressed 
and the chiral properties of the Dirac operators used in this thesis are out- 
lined. As indicated above we define the theory on a lattice with boundaries. 
The specific form of the boundary conditions and the implications of the 
presence of the boundaries on the renormalisability are discussed thoroughly 
in Chapter [H The two dimensional fermionic model we utilise in this work 
is carefully defined in Chapter O In particular, we are concerned with the 
symmetries of the model, its lattice formulation and renormalisability. Since 
Wilson fermions explicitly break chiral symmetry, it has to be assured that 
the symmetry is recovered in the continuum limit. The employed strategy 
and the result are presented in Chapter [6l We use this result in Chapter [7] to 
define renormalised couplings. A next-to-leading order computation is then 
carried out for Wilson and Ginsparg- Wilson fermions. We draw conclusions 
and give an outlook in Chapter [3 



Chapter 2 

Lattice perturbation theory 



In this Chapter we introduce the basic concepts used in this thesis. Since 
we want to use lattice perturbation theory, we have to discuss perturbative 
renormalisation (Section I2.1I) . After these general remarks we list our no- 
tation and conventions for the lattice computation (Section 12.21) . We close 
the Chapter with some remarks on the continuum limit and the analysis of 
lattice diagrams (Section l2.3p . 

2.1 Renormalisation 

The bare couplings and masses that appear as parameters in the classical 
action of a quantum field theory are not the couplings and masses which 
are measured in experiments. Experimentalists rather gather data of cross 
sections and transition amplitudes. These quantities have to be computed in 
the theory that is supposed to describe the phenomenons. For each parameter 
in the action one input measurement is needed. Once all the parameters are 
fixed the theory can be used for predictions. 

The crucial point of course is how many parameters are there in the 
action. The more parameter the less predicting power does the theory have. 
The number of terms in an action and thus the number of bare parameters 
is mainly restricted by symmetries and dimensional analysis. 

Computing a cross section or transition amplitude yields a relation be- 
tween an observable and the bare parameters of the theory. The observable 
itself may now be called coupling. In order to avoid confusion one calls it 
renormalised coupling since it is a redefinition of the bare coupling. In the 
same way all other couplings and masses may be redefined. The renormalised 
quantities may be regarded as the physical parameters of the theory because 
all observables can be expressed in terms of them. 
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In perturbation theory the necessity for renormalisation is encountered in 
the form of ultraviolet infinities when calculating loop corrections. To keep 
physical amplitudes finite these infinities have to be absorbed in a redefinition 
of the parameters and fields order by order in the perturbative expansion. 

In order to handle the terms producing the infinities they first have to be 
rendered finite. In lattice perturbation theory the inverse lattice spacing 1/a 
provides an ultraviolet cut-off to the theory. This regularisation has to be 
removed before comparing with experiment. On the lattice this amounts to 
taking the continuum limit a ^ 0. In this process the ultraviolet divergences 
show up and the renormalisation has to be implemented. 



2.1.1 Mass independent renormalisation scheme 

Let us consider a quantum field theory with one mass and one coupling con- 
stant that has been regularised on an infinite lattice, say QCD with N mass 
degenerated quarks. All information of the theory is contained in the n-point 
Green's functions. Any unrenormalised Green's function T{p; go, mo, 1/a) will 
then depend on the momenta of the external lines collectively labelled with p, 
on the bare mass mo and coupling constant go and on the ultraviolet cut-off 
1/a. 

QCD is a renormalisable quantum field theory. A renormalisation scheme 
is given through conditions that define the renormalised mass m/j, renor- 
malised coupling gu. Often also a wave function renormalisation factor Zi 
for each type i of fields in the theory (i.e. in QCD one for the quark fields and 
one for the gluon fields) is introduced. This is not a necessity but convenient 
in the course of computations. Since we will consider massless field theory, 
with a renormalisation scale /i, a mass independent renormalisation scheme 



is needed to avoid infrared divergences [35||. The renormalisation conditions 
are then posed at the scale fj, and vanishing renormalised mass. For the 
renormalised parameters one expects 

Zi = Zi{go, afi) , (2.1) 
gR = goZg{go,afi) , (2.2) 
niR = nig Zm{go, an) , mg = mo-mc, (2.3) 

where rric accounts for the additive mass renormalisation needed if chiral 
symmetry is broken by the regularisation (cf. Section 13. 4p . If the regulari- 
sation does not violate chiral invariance the renormalised mass vanishes at 
zero bare mass. 

The renormalized coupling in lattice QCD, for example, may be defined 
through demanding the triple gluon vertex function to take its tree level value 



2.1 Renormalisation 
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at momenta of order /x. Then the renormalised Green's functions have finite 
continuum limits. They are functions of the renormalised coupling and mass 
and are related to the bare ones as 

TR{p;gR,mR,ix) = Zr{go, aij)T{p; go,mo,l/a) , (2.4) 

where Zr depends on the number and types of the external lines. (Note that 
there will also be some dependence on a gauge fixing parameter as in the 
continuum. However, this dependence disappears when physical quantities 



are computed and is not important for the aspects considered here. See |36|| 
for a complete review of lattice perturbation theory.) 

Eq. (|2.4I) really only holds in the continuum limit. At finite cut-off, that 
is finite lattice spacing a, perturbation theory states that the renormalised 
Green's functions are cut-off independent only up to terms of order a 

Tr{p; gR, niR, /i, afi) = Tr{p; gR, tur, jJ,) + O (a(ln a)'') , (2.5) 



at fc-loop order [6||. These terms are called scaling violations. Since they are 
small near the continuum limit we suspend their discussion until Section 12.31 
and neglect them in the following. 

2.1.2 Renormalisation group equations 

Since gR and rriR depend on the renormalisation scale /i while T does not, 
differentiation on both sides of (12. 4p yields the so called renormalisation 
group equations 

where 

d 

PigR) = f^-gj^gnigo, a^) , (2.7) 
d 

T{gR) = In ^m(fi'o, a-lA , (2.8) 
d 

IrigR) = fJ'-gj^ In Zr{go, a/i) . (2.9) 

The renormalisation group functions /3, r and 7r are the so called beta- 
function for the coupling and the anomalous dimensions of the mass and the 
Green's function. (If two types of fields appear in F, say rii of type one and 
n2 of type two, we have to take Z^ = Z'l^Zl^^ . Then 7r is the sum of the 
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anomalous dimension of the two types 7r = ^171 + ^272-) Note that the 
coefficient functions f l2.7H2.9l) must be independent of a because they appear 
in a differential equation of an cut-off independent quantity. Since they are 
dimensionless they must also be independent of /i. Thus they only depend 
on the renormalised coupling gj^. 

The functions (3, r and 7^ can be calculated in perturbation theory as 
a power series in the renormalised coupling. In the case of QCD the beta- 
function 

f3{as) = b^a] + ha] + 62^^ + O(a^) , (2.10) 



is known to tree loops. The first two coefficients, for example, are [37|, 13; 



39,^ 



,„._(n-f) ,..-(l02-5|^^). ,2.11) 

(Note that in many textbooks and publications another definition of the 
renormalised coupling of QCD is used. The relation to the one used here is 
«s = ^y47r.) 



2.1.3 Asymptotic freedom 

From the shape of the beta-function the behaviour of the renormalised cou- 
pling at high energies may be deduced. The example above is characterised 
by a negative [3{g) for small 5^ > and leads to a vanishing renormalised 
coupling as /i — ^ cxD. This behaviour is called asymptotic freedom. There 
are three other possible scenarios, we do not list them here but refer to the 



diverse textbooks on the topic [4ll, l42, l43|. To make the above statement 



more explicit and general, assume a beta-function that is negative for small 
positive g 

m'^'-b9\ b>0, (2.12) 

where g"' is the power of the coupling in front of the lowest-order divergent 
diagram contributing to P{g) and therefore is always greater one. The renor- 
malisation group equation is then 

/i^^?(/") = -bg'^ifi) . (2.13) 

Given the renormalised coupling at some scale /i the coupling at the energy 
E can be calculated by integrating this equation. The solution is 

giE) = gi^^) [1 + (n - 1) 6 ln(E//i) (gifi))^-"'^] "'^^""'^ . (2.14) 



2.1 Renormalisation 
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For E ^ oo this solution becomes independent of g{fi) 

g{^E) [{n - 1) h ln(E/^)]^'/("'') . (2.15) 

Thus starting from a value justifying the approximation (I2.12p g{E) always 
tends to zero for ^ oo. On the other hand at small E the coupling may 
become large g{E) > 1. Thus perturbation theory becomes unreliably at 
small energies and non-perturbative methods are needed. 

Along similar steps it can be shown that the effective dimensionality 
of operators and fields is given by dimensional analysis up to logarithmic 
corrections 
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In this context it is worth mentioning that the first two coefficients of 
the beta-function are independent of how exactly the renormalised coupling 
is defined as long as for small bare coupling Qr = Qq + 0((y'o)Q To see this 
assume two renormalised couplings qa and qb- Since there are no other 
dimensionless parameters qa is a function only of qb- We can expand the 
one in powers of the other 

gAigs) =gB + c^gl + 0{gl) , (2.16) 

or 

gsigA) =gA- cig\ + 0(^7^) , (2.17) 

where the leading order coefficient is fixed by the condition that qa and qb 
at leading order are equal to the bare coupling. The two beta-functions can 
be related through 

, d esSI dgB dgA m n f ^ 9gA ^ 

PAigA) = f^^^gA = /.^^ = PBigB)g^^. (2.18) 

The beta-function has an expansion in the renormalised coupling 

PsigB) = b^gl + bfgl + 0{g%) , (2.19) 

where the leading power is two but the argument holds for arbitrary leading 
power greater unity. In terms of gA this becomes 

/3b(^?a) = b^gl + (b? - 2c,bl^)g\ + 0{g\) , (2.20) 
and the derivative is 
dgA 



dgi 



1 + 2c^gB + 0{g'B) = 1 + 2c^gA + 0{gi) . (2.21) 



^Note that notation might be misleading here. With a renomaUsed coupUng like as 
of QCD is meant. The more familiar renormlised coupling of QCD, that is called ga, the 
beta-function woud start with a third power. See also the note under (|2.1ip 
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Now the right hand side of (12.18^ can be evaluated in terms of qa 

Pa{9a) = [bEgl + (bf - 2c,b^)9\ + 0{g\)] ■ [l + 2c,gA + 0{g 



(2.22) 
(2.23) 



proving that the first two coefficients are universal in the sense that they 
neither depend on the regularisation nor on the the renormalisation scheme. 
For the first coefficient this is a direct consequence of demanding the renor- 
malised couplings to coincide at leading order (that is, at leading order they 
coincide with the unrenormalised coupling). The second order coefficients 
coincide because the expansion of beta-function starts at the next-to-leading 
order. 

Finally we introduce the A-parameter 



A = f^ibogl) 



-6i/(262) -l/(26o4) . 



exp 



dg 



Pig) 



+ 



bog^ big 



(2.24) 

In the massless theory the A-parameter is the only dimensionfull parameter. 
It is the standard solution of the renormalisation group equation for physical 
quantities 

/^^ + /5(^7ii)^j P{fi,gR)=0, (2.25) 

which expresses the abitrariness of the refernce scale fi. 

For each the two renormalised couplings gA and gs A-parameter can 
be defined. Given the relation (12.161) between the couplings the ratio of the 
A-parameters is then a pure number 
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Aa/A 



B 



exp 



2^0 Wb 9a 



(2.26) 



2.1.4 Multiple couplings 

So far we considered theories with a single dimensionless coupling. It is not 
difficult to generalize the concepts to multiple such couplings |4ll|- There will 
be as many renormalised couplings gi as bare couplings. Green's functions 
depend on all these couplings and in (|2.4I) gi^ and g^^ may collectively refer 
to them. Then for each gi there is a renormalisation group equation 



A(^(/^)) , 



(2.27) 



with /9; depending in general on all the renormalised couplings gi. In the case 
of one coupling the beta-function determines the asymptotic behaviour of this 
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coupling. In the case of multiple couplings the beta-functions /?/ determine 
the asymptotic trajectories in the space spanned by the couplings gi. Clearly, 
there are many possibilities now. Let us concentrate on the prominent case 
of trajectories approaching a fixed point in (7-space. A fixed point = 
is defined through a mutual zero of the beta-functions 

Piig*) = 0. (2.28) 

Shifting gi ^ gi — gi by the fixed point, Taylor-expanding Pi{g — g*) and 
ignoring terms 0{{g — g*Y) eq. (12.271) becomes 

/^^[^?K/i)-^7r] =5^ Mfc[(7,(/i) (2.29) 

^ k 

with the matrix M given by 

M„ = ( ^) . (2.30) 



dgk 



9=g 



Suppose that the eigenvalues of this matrix are non-degenerate. Surely, that 
is not always true but it is the generic case. Then the eigenvectors v"^ 

5^Mrf = A>'", (2.31) 

k 

form a complete set and can be used to express the solution to (12.29P 

giifi) = g: + J2c^vrfi^-, (2.32) 

m 

with coefficients Cm- 

The qualitative behaviour for /i ^ cxd is thus governed by the eigenvalues 
A™ and the coefficients Cm- In particular, the fixed point is approached if 
and only if = for all A™ > 0. A zero eigenvalue may be caused by a 
vanishing in which case the fixed point is reached for any value of this 
coupling in a region around the fixed point. 

The eigenvectors of the negative eigenvalues define a subspace containing 
the trajectories attracted by the fixed point. Trajectories with support out- 
side of this subspace may get very close to the fixed point but are eventually 
repelled. 

We close this section by pointing out that the eigenvalues A™ are in- 
variant under a change of basis, that is going to a differently defined set of 
renormalised couplings gi. The change g ^ g amounts to a similarity trans- 
formation of the matrix M and thus preserves its spectrum. This can be 
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seen as follows. The new couplings will be functions of the gs and satisfy 
renormalisation group equations 

f^-^Uf^) = E = Piim) ■ (2.33) 

Thus (3 transforms as a contravariant vector in coupling space 

m-j:^^rn{9). (2.34) 

Differentiating on both sides with respect to we get 

E§^%^ = E/f-«9) + Et^— ■ (2-35) 

^ ogm ogk ^ ogmogk ^ dgm gk 

At the fixed point g* the first term on the right hand side vanishes and we 
are left with the matrix equation 

MS^SM, (2.36) 

with 

^ ^ 9=9(9 ) ^ ^ 

As long as 5" is invertible this is a similarity transformation and the eigen- 
values of M are those of M. 



2.2 Discretisation 

The regulator used in the computations of this thesis is the Euclidean lat- 
tice. If the (Minkowskian) time coordinate is Wick rotated to imaginary 
(Euclidean) time 

= ix^ , (2.38) 

the imaginary unit in front of the Minkowski-space action in the path integral 
of a quantum field theory becomes a minus sign 

J D e^^^ ^ / ^ ^^^'^ ■ ^^-^^^ 

If the action is bounded from below (which is the case for physically relevant 
theories) the weight factor can be interpreted as a probability distribution for 
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field configurations. Evidently there is a close connection between field theory 
and statistical physics if the weight is identified with the Boltzmann factor. 
The only subtle point at this stage is, one has to assure that the analytic 
continuation of the n-point Green's functions to imaginary time exists (see 



Section 1.3 in [42| and references therein). The path integral becomes a 
mathematically well defined object, that is a convergent multidimensional 
integral, by discretising Euclidean space-time on a finite lattice. This is the 
basis of non-perturbative Monte-Carlo techniques. 

From now on we work in D = d+1 dimensional space-time with Euclidean 
metric and drop the superscript, that is Xq refers to Euclidean time (Greek 
subscripts always run from to d). The Euclidean Dirac matrices satisfy 
anti-commutation relations 

{7^, 7,} = 25^,, (2.40) 

are all hermitian 

7i = 7m , (2.41) 
and are related to their Minkowskian counterparts as 

7o = 7o^', 7^ = -^l■'■ (2.42) 

We are here dealing with theories in two and four dimensions and the defi- 
nition of the Euclidean 75 differs by a factor due to demanding hermiticity 

D = 2: 75 = z7o7i , D = 4 : 75 = 7o7i7273 • (2.43) 

Explicit representations of the Dirac matrices are given in the Appendix 

Although the ultimate goal is to consider field theories on a lattice with 
boundaries we start here with the more common hypercubic lattice with 
either infinite extension or periodic boundary conditions. The sites of the 
lattice are labelled by = an^ with integer and lattice spacing a, which 
is the same in all directions. On a finite lattice the coordinates are restricted 
to < < L giving a total number of lattice sites V/a^ = /a^ . 

Continuum space-time integrals are on the lattice replaced by sums over 
all lattice sites 

' d^x ^ a^'J^- (2.44) 

X 

The lattice spacing a is the minimal distance in the system and thus in- 
troduces an ultraviolet cut-off. The momenta can be restricted to the first 
Brillouin zone 

- - < < - • (2-45) 
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On an infinite lattice continuum momentum integrals are cut-off 



On a finite lattice the allowed momenta are a discrete set in the range (I2.45p . 
For periodic boundary conditions and integer the V/ allowed momenta 
are 

2nn., , , , , 

Pt. = , = -L/2 <n^< L/2 , (2.47) 

and the momentum integrals also become momentum sums 

77 E- (2-48) 



7r/a 



(27r)^ V 



p 



The discretisation of integrals was straightforward. However, continuum 
differential operators have infinitely many valid lattice representations. We 
introduce here the simplest possibilities which will serve as building blocks 
for more difficult choices. We consider lattice fields ip{x) defined at the sites 
of the lattice. On a finite lattice we have to specify boundary conditions. We 
choose general periodic boundary conditions 

il){x + Ljl) = e'^^^^ix) , -n <e^<n, (2.49) 

parametrised by the phases 6^^. The forward and backward finite difference 
operators are 

d^ij{x) = l[ij{x + afi) - ij{x)] , (2.50) 
^;^Pix) = i[V'(x) - ^(x - afi)] , (2.51) 

where /i is a unit vector in /^-direction. 

There is a different way of incorporating such general boundary condi- 
tions. One takes the lattice fields as periodic 

i){x + Lp.) = ip{x) , (2.52) 

and defines the forward and backward finite difference operators as 

d^ipix) = ^[Xf^ipix + afi) - ip{x)] , (2.53) 
dlij{x) = \[i,{x) - \-^i,{x - afi)] . (2.54) 

The phase factors depend on 

= e^"^"/^ , -TT < < TT . (2.55) 



2.3 Continuum limit and lattice artefacts 
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The two notations are connected by an Abelian gauge transformation. A 
non-zero introduces a "momentum" that is not restricted to the values of 
(12.47p . The construction with the phase factors in the finite differences is 
computationally easier and therefore adopted here. On infinite and periodic 
lattices the forward and backward finite difference operators obey 

id^)^ = -d;. (2.56) 

Other important lattice operators are the anti-hermitian averaged finite 
difference operator 

d^ilj{x) = + d;)ij{x) = i^[XM^ + afi) - X;'^{x - aft)] , (2.57) 

and the hermitian lattice Laplace operator 

d^dli){x) = d^dli){x) = j^iXfiipix + afi) + \~^^ix - aft) - 2ijj{x)] . (2.58) 

As already mentioned there is some freedom in discretising differential 
operators and therefore the action of a given field theory. This is due to 
the smaller symmetry on the lattice. In particular Lorentz invariance is 
broken and an infinite number of (in this case irrelevant from the point of 
view of renormalisation) terms can appear in the lattice action. In this way 
infinitely many different lattice actions for the same continuum theory are 
possible. However, they are expected to be equal in the continuum limit, 
that is when the cut-off is removed. One says the different lattice actions 
fall into the same universality class characterised by the target continuum 
theory. 

In numerical Monte-Carlo computations one is interested in lattice actions 
that balance between complexity (numerical cost) and the rate at which 
the continuum limit is approached (systematic error). Lattice perturbation 
theory is an essential tool to provide analytic understanding of these lattice 
artefacts. 



2.3 Continuum limit and lattice artefacts 

The renormalisation group function P{gR) describes the variation of the 
renormalised coupling gji with the cut-off at fixed bare coupling 



Pig. 



R) 



/i 



dgRigo,aiJ^) 
djji 



90 



da 11 



(2.59) 



90 
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Now we can ask how the bare coupling has to be varied with the cut-off for 
fixed renormalised coupling. The lattice beta- function /3lat(5'o) is defined 
through 

d ( d d ] 

"da^^ " [da ^latIa-o)^ | 5'i?(fl'o, a/i) = 0. (2.60) 



Plat{9o) = -a- 



da 



f2.611 



Since 



d d 

a—gR{go, afi) = fi—gn^go, afi) = (3{gR) , (2.62) 

the two beta-function are related in the following way 

/5lat(ao) ^ = f^ign) ■ (2.63) 

ago 

And because at lowest order in perturbation theory the both couplings are 
equal, gn = go + O^g^), we know from Section 12.1.31 that the first two 
coefficients of the beta-functions are identical. In particular, this means that 
the bare coupling vanishes in the limit a — > 0. 

As mentioned at the end of Section 12.1.11 renormalised lattice Green's 
functions have a finite continuum limit a —>■ and differ from this limit by 
terms of order a 

Tr{p; gR, niR, /i, a/i) = Tr{p; gR, tur, /i) + O (a(ln a)'') , (2.64) 

at /c-loop order in perturbation theory It is widely believed that this be- 
haviour also holds beyond perturbation theory and non-perturbative Monte 
Carlo data (naturally only available at finite a) is extrapolated to the con- 
tinuum accordingly. Since perturbative computations are the only analytic 
tool to learn about the size of these lattice artefacts such computations are 
essential, especially if new methods are used. For example, the amplitude of 
the lattice artefacts can be very different for different lattice actions. 

In order to obtain numbers in lattice perturbation theory one almost 
always is forced to evaluate some lattice or momentum sums numerically. 
Consider a quantity P{a/L) that is a sum of lattice diagrams at 2-loop order, 
where we suppress any dependence on the external lines. We supose that P 
is dimensionless. If it is not, it can be made so by appropriate factors of 
a. Futhermore we suppose that it is finite at a/L = 0, which always can be 
achieved by multiplication with appropriate factors of a/L. Then P has an 
expansion in a/L 

oo 

P{a/L) = Y, [rn + Sn Ha/L) + t„ \n\a/L)] {a/LY . (2.65) 

n=0 
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The generalisation to abitrary loop order should be obvious. 

In practice, P{a/L) is computed at several values of a/L. In order to 
extract the coefficients of the expansion 12.651 we use the method described 
in Appendix D of Ref. [i^. In this way it is possible to reliably determine 
the systematic uncertainties that are inevitable involved in such a numerical 
analysis of lattice Feynman diagrams. 
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Chapter 3 

Chiral symmetry on the lattice 



The Euclidean Lagrangian of free fermions is 

£ = ^(x) (7^(9^ + m)V^(x) , (3.1) 

where the fermion fields carry suppressed Dirac and flavour indices. The 
flavour indices, labelling the N fermions, are contracted by a unit matrix in 
flavour space. For m = the theory is invariant under a global U (N) x U{N) 
flavour symmetry which can be decomposed into U{l)v x U{1)a transforma- 
tions 

acting equally on all flavours and chiral SU{N) x SU (N) transformations 

where the generators of the SU (N) algebra A"^ act on the flavour indices. 
The subscripts V and A refer to the associated Noether currents (see next 
section), which are of Lorentz vector and axial- vector type: 

= ^-fi^tp , = V'T/.TsV' , (3.4) 

and 

\/; = V^7/.A>, A^ = V^7m75A>. (3.5) 

We use the same letters for the SU (N) singlet (13.41) and vector currents (13. 5p 
but indicate the difference by an additional superscript for the vector ones. 
We can introduce a new set of generators 

^2 = (l-75)A^ t^ = (l + 75)A% (3.6) 

^The generators Xa are normalised to obey Tr {A° A*"} — 2Sab- See Appendix I A . 1 .21 for 
more details. 
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with commutation relations 

[tl,ti]=2zr'Hl, (3.7) 

r«,4]=2zr%, (3.8) 

[tlA]=0. (3.9) 

The tis and tns obviously form two closed subalgebras. Thus, as indicated 
by the notation, chiral SU (N) x SU (N) is the direct sum of two SU (N) 
subgroups acting independently on the left- and right-handed components of 
the fermion field 

V^L = |(1-75)V^, ^pR = 1(1 + 1,)^. (3.10) 

A general mass term ipmip with m = diag(mi, ...,771^) in (13.11) breaks 
all of these symmetries but multiplication with an U{l)v phase. N mass 
degenerated fermions (m oc 1) lift this to an U{N)y since also the generators 
A" remain unbroken. 

In QCD chiral symmetry plays a key role in understanding the mass 
spectrum of the light mesons. The pious, for example, are seen as the Gold- 



stone bosons [46 
SU{2) X SU{2) 



47l | associated with the spontaneously axial generators of 
48|. The QCD Lagrangian with only the two light u and d 
quarks has this symmetry in the massless limit, which is a good approxima- 
tion because m„ ^ is much smaller than a typical hadron mass scale. Indeed 
isospin symmetry and quark number conservation, U{l)v and SU{N)v re- 
spectively, are experimentally confirmed to high precision [l|. The axial 
generators 75 A"^ are spontaneously broken in the quantum theory by a non- 
vanishing quark condensate leading to three massless Goldstone bosons: the 
pious. The small but non-zero mass of the pious is due to the fact that 
SU (2) X SU{2) is only an approximate symmetry. 

So far, nothing has been said about the axial U{1)a^ i.e. continuous 
chiral phase transformations ip —>■ e^'^^'^^ip. In QCD this symmetry is also 
broken, but it must be in a different way since there is no light iso-singlet 
state expected from the Goldstone theorem. The effect can be traced back to 
the topological structure of the vacuum gauge field ji^ and also explains the 
suppression of the electromagnetic decay of the neutral pion [i^. However, 
this chiral or axial anomaly is tightly connected to gauge symmetry. The two- 
dimensional fermion model considered in this thesis has no gauge symmetry 
and therefore the U{1)a is not anomalous. 

As already stated above understanding the structure of the low energy 
regime of QCD needs understanding of chiral symmetry and how it is vio- 
lated. Since the QCD coupling is large (cf. Section I2.1.3P at low energies 
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non-perturbative methods such as numerical lattice QCD are needed. The 
discretisations used today handle chiral symmetry quite differently. We will 
discuss two different lattice Dirac operators and present their chiral prop- 
erties. But before that we introduce the currents associated with the sym- 
metries and derive operator identities from infinitesimal variable transforma- 
tions in the path integral. 



3.1 Continuum Ward identities 

Here we derive operator identities, so called Ward identities. They are con- 
veniently derived by variable transformations in the path integral for the 
expectation value of operators, where the new variables are connected to the 
old ones by infinitesimal local transformations 

i^ip + di, i^i + d^. (3.11) 

The expectation value of a local operator C, that is a field composed of 
fermion fields and their derivatives evaluated at the same space-time point, 
is 



(0> = i j D^D^Oe-^ (3.12) 



with the Euclidean action 

]0™ r — I aD 



S = yd^x£ = y d^x^(x)(7^(9^ + m)^(x). (3.13) 

A linear change of the Grassmann valued fermion fields in f l3.12p not 
only effects the operator and the action but also the integration measure by 
introducing a non-trivial Jacobian 

ip' = Aip, i' = Ai, (3.14) 

B^Bip = det AdetA- B^'Bip' = J ■ Bi'Dip' . (3.15) 

For simplicity we assume the infinitesimal transformations (I3.1ip can be rep- 
resented by 

A = 1 + ujX, A = 1 + ujX, (3.16) 

where u is an infinitesimal function of space-time evaluated at the point of 
the field and X , X are matrices acting on Dirac and flavour indices. Then 

we can use det(l + u;X) = 1 + ujTt {X} + 0{uj^) to write J = l + 6J + 0{uj'^). 
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Let 6S and 60 be the change in the action and the operator for a given 
change of fermion fields. Performing an infinitesimal change of variables in 



the path integral yields 

(0> = | J B^D^Oe-', (3.17) 

= ^ jBiiDij{l + 5J + 0{uj^))iO + 5O)e-^{l-5S), (3.18) 

= (O) + {60) - {06S) + {05J) + 0(cj2) . (3.19) 
Thus we derived the general identity 

{50) - {05S) + (C 5J) = . (3.20) 



Let us restrict uj{x) to a region TZ, that is uj{x) = for a; ^ 7^. If we choose an 
operator Oext defined outside the region IZ where the variable transformation 
is performed then 5Cext = 0. Likewise the change in the action is supported 
only in IZ 

SS= [ d^x6C. (3.21) 
Jn 

Eq. f l3.20p further simplifies to 

(axt55> = (axt5j>. (3.22) 

Now we promote the global transformations of (j3.2l) and (13. 3p to local 
but infinitesimal transformations like 



a) 


u{l)v: 


Sip = 


iUylp , 




—iuvi' 


(3.23) 


b) 


u{1)a: 


dip = 


iujAl^i^ , 




iujAi^lb 


(3.24) 


c) 


su{N)v: 


6ip = 


tUJ^X^lP , 




-iu^^X'' 


(3.25) 


d) 


suiN)A: 


5ip = 








(3.26) 



where the us are now infinitesimal functions of space-time evaluated at the 
point of the field. 

For the SU{N) vector transformations c) and d) the Jacobian is unity 
{6 J = 0) because the SU{N) generators are traceless. The change in the 
action is 

SvS= [ d^xu;^[-9^V;"(x)-V^(x)[A",m]V^(x)] , (3.27) 
Jn 

and (3.28) 

S^S = [ d^x [-d^^A^ix) + 2mP'^(x)] , (3.29) 
Jn 
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respectively. Then the identity (13.221) becomes 

{d,V;{x) Oext) = (V^(x)[m, X-],lj{x) Oext) (3.30) 

and 

{d,Al{x) Oext) = (?(x)75{A^ m}iP{x) Oext) • (3.31) 

In deriving these identities one has to make use of the arbitrariness of the 
exact definition of a;(x). For mass degenerate fermions one finds the SU{N) 
vector current conservation 

(a^F;(a;)axt> =0 (3.32) 

and the partially conserved axial current (PCAC) 

(9^A;:(x) Oext) = 2m{P%x) Oext) , (3.33) 

with the flavour vector pseudo-scalar density 

P"(x) = V^(a;)75AXx) . (3.34) 

Although we considered here free fermions witho ut g auge interactions the 
calculation in QCD yields the same identities (see [49j for example). 
For the U{1) transformations a) and b) the change in the action is 

6vS = I d^x uoy [-d^V^{x)] , (3.35) 
Jn 

and (3.36) 

5aS = / d^x uja l-d^A^ix) + 2mP{x)] . (3.37) 

The corresponding operator identities are the fermion number conservation 

(9^V;(x)axt> = (3.38) 

and the singlet PCAC relation 

(d^A^ix) Oext) = 2m{P{x) Oext) , (3.39) 

with the singlet pseudo-scalar density 

P{x) = i^{x)-f5ij{x) . (3.40) 

As already mentioned in QCD the U{1)a symmetry is anomalous and 
Afj_{x) is not conserved in the quantum theory even for vanishing masses. 
Along the presented steps, also in QCD we would have arrived at (13.391) . in 
contradiction to the anomaly. This is because the treatment here was very 
formal. That is, we performed variable transformations in an integral that is 
not well defined in first place. If QCD is regularised first (13.391) is changed 
and the anomaly is recovered. Using Ginsparg- Wilson fermions the anomaly 
can be computed easily [50| . 
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3.2 Lattice Ward identities 

The steps that led us to the continuum PCAC relations can be applied on 
the lattice [Slfl- The expressions for lattice currents are then multilocal, 
in the sense that they involve more then one lattice site, and additional 
terms can appear that are allowed by the lattice symmetries. However, the 
difference between multilocal and loca@ becomes irrelevant in the continuum 
limit. Futhermore the continuum PCAC relations are part of the definition 
of the theory and thus, at finite lattice spacing and for properly renormalised 
operators they hold up to cutoff effects. 

In perturbation theory on expects terms of 0(aln(a)') in /-loop order. 
In the case of massless fermions, for example, eq. (13.331) implies the lattice 
version 

{d, {A;{x))^ (axt)R> = O(aln(a)0 . (3.41) 

On the lattice we also use the local currents and densities. Thus the contin- 
uum definitions (EIH]), ((331), (EIHll), flXiOll transfer to lattice. Although the 
axial-vector and vector currents are conserved in the massless theory they 
receive a finite renormalistion on the lattice 

{A;{x))^ = ZaA';{x) and iA,{x))^ = ZaA,{x) , (3.42) 
{V;ix))^ = ZyV^ix) and {V.ix))^ = ZyV.ix) , (3.43) 



with 



Zi = l + Zf" g^ + ... . (3.44) 



3.3 Nielsen-Ninomiya theorem 

Before discussing the two lattice Dirac operators used in this thesis we briefly 
describe the broader context of lattice fermions and chiral symmetry. As we 
will see chiral symmetry in the lattice regularisation is tightly connected to 
other desirable properties of the operator D = ''i^D^j, in the massless action 

S = a^^^{x)Dij{x). (3.45) 

X 

^Here locality refers to operators composed of fields that are all taken at the same 
space-time point. It may not be confused with the locality of a differential operator like 
the Dirac operator. 
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The kernel D{x — y) and the Fourier transform D{p) of the Dirac operator 
are defined through 

Dij{x) = 5^ D{x - y)ij{y) , D{x - y) = T^" e^^(^-^)5(p) . 

(3.46) 

To give an example consider the simplest lattice Dirac operator one could 
think of 

Dn.i.e = lbM + d^)}^ (3-47) 

that is the averaged finite difference (12.571) contracted with the gamma- 
matrices. Indeed, -Dnaive anti-commutes with 75 and thus (I3.45p would be 
invariant under the chiral transformations (13.241) and ( 13.26^ . But it turns 
out that this Dirac operator leads to a proliferation of fermion species. The 
operator in Fourier space 

5naive(p) = Z XI ^^"^(^Pf^) ' ^^/^ <Pti<'^/a', (3-48) 



has 4 and 16 zeros in two and four dimensions respectively. Thus D~^^^^ is 
not only propagating a single physical fermion (the zero at = 0) , but also 
the so called doubler modes. 

This is not an accident. The appearance of the doublers is in accordance 



with the Nielsen-Ninomiya theorem [52, 153|, 15J, 155|. It can be stated in the 
formic 

Theorem 1 Any massless lattice Dirac operator cannot satisfy the following 
four properties simultaneously: 

a) aD(x) is local in the sense that it is bounded by (7e~'^'^^^'/'^ 

b) D{p) = i'j^Pf, + 0{ap^) for \p\ <^7i/a 

c) D{p) is invertible at all non-zero momenta 

d) 75^ -D^, = 

In the bound in a) C and 7 > are constants that do not depend on a. 
With respect to the universality of the continuum limit a Dirac operator 
with such exponentially small tails is certainly as good as an ultra-local 
operator with only nearest neighbour interactions. The property b) is the 



^The theorem was originally proven for fermions on two- and four-dimensional lattices 
0. 
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right continuum behaviour at small momentum. Doublers are excluded by 
c) and chiral symmetry in the form of ( 13.24^ and (I3.26p is guaranteed by d). 

Any fermion discretisation in two or four dimensions has to abandon at 
least one of these desirable properties. In the example given in this section 
one finds a bunch of doubler modes, i.e. c) is violated. In the free theory this 
may not be a problem, but in the interacting theory they would contribute 
in the loop corrections. 

In the following two sections we discuss the Wilson and the Ginsparg- 
Wilson fermions, both abandoning d), but in a very different way and with 
very different consequences. 



3.4 Wilson fermions 

In the view of the last section Wilson's approach (s^ to remove the doublers 
is to sacrifice chiral symmetry by adding an irrelevant term to (13.471) . The 
Wilson-Dirac operator is given b}0 

= \{l^.{dl + d^) - ad*M , (3.49) 

and corresponding massive Dirac operator is simply 

Dm = -Dw + mo . (3.50) 

Looking at the Fourier transform 

^w(p) = i 5Z{7Msin(ap^) + f sm\ap^/2)} , (3.51) 

we see that the doublers receive a mass of the order of the cut-off and hence 
are strongly suppressed at finite a and eventually disappear in the continuum 
limit. 

The price to pay is that the lattice Laplace operator adj^d^ explicitly 
breaks chiral symmetry at finite lattice spacing. As a consequence a vanishing 
bare mass does not automatically imply a conserved axial current as in the 
continuum (13.330 . But this operator identity has to hold for the renormalised 
quantities also on the lattice up to scaling violations 

(axt)R> = O(aln(a)0. (3.52) 



''Often there is an additional parameter < r < 1 multiplying the lattice Laplace 
operator: —ard*^d^. Throughout this thesis we set r = 1. 



( d,{Al{x))^-2mn{P{x)) 
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In fact, together with (12.30 

tur = ruq Zraigo, an) , mg = mo-mc, (3.53) 

eq. (I3.52P serves as definition of the critical mass nic [51], that is, the value of 
the bare mass at which the axial current is conserved up to scaling violations 

mo = mc{g) such that (d^ (^m(^))r (^ext)R> = 0(aln(a)') . (3.54) 

In the two-dimensional Gross-Neveu model the global chiral U{1) symmetry 
is not anomalous and thus a critical mass for the conservation of the singlet 
axial current can be defined 

mo = TUcig) such that {df, {Af,{x))^ (Cext)R) = 0(aln(a)') . (3.55) 

In this thesis we employ the latter definition. As indicated the critical mass 
can be computed in perturbation theory as a power series in the coupling 
constant 

Having mass dimension one, rric will cancel the linear divergence spoiling 
the PC AC relation (( 13.33p or (I3.39p ) due to the explicitly broken chiral sym- 
metry. However, there may also be a term of order 0(1). In QCD such a 
term is not present because there is no chiral symmetry breaking operator of 
mass dimension four. Therefore it is enough to tune m^. But in two dimen- 
sions there are chiral symmetry breaking four fermion interactions and one 
has to tune a dimensionless parameter in addition to nic (see Section [5.3. II) . 



3.5 Ginsparg- Wilson fermions 



From the Nielsen-Ninomiya theorem (page 1251) one may conclude that it 
is not possible to construct a meaningful lattice theory of fermions with 
chiral symmetry. However, the way out was discovered shortly after the 
original paper by Nielsen and Ninomiya. Studying block-spin renormalisation 
in lattice QCD, Ginsparg and Wilson 33 1 found the relatioE@ 



j,D + D^, = aD^,D, (3.57) 

for the lattice Dirac operator. But it was not recognised at the time and 
remained unappreciated until Hasenfratz [57| "rediscovered" it in the context 



®In Ref. [33| the authors actually derive a more general relation, but (|3.57p is the most 
spread version and is used solely in this thesis. 
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of the perfect action. Indeed, the Dirac operator of the perfect action obeys 
the Ginsparg- Wilson relation (13.5711 . 

In the following Neuberger [58|| found an explicit solution to ( 13.57p and the 
dimensional reduced Dirac operator of the domain wall fermion was recog- 



nised to obey the Ginsparg- Wilson relation [59 



Before turning to the solution presented by Neuberger we point out that 
the importance of (13 .57^ is due to the fact that it allows for an exact chi- 



ral symmetry |50l |. Consider the infinitesimal transformation of the fermion 
fieldB 

^ ^ + e^^il - aD)^ , V'^V + e^Ts- (3.58) 
This is an exact symmetry of the action (13.451) 

55 = ^ e^{x) (^75(1 - aD) + 75/^) ^(x) ^ . (3.59) 

X 

In the continuum limit (I3.58P become the familiar chiral phase transforma- 
tions (I3.24p . Chiral flavour transformations are deflned analogously by re- 
placing e75 e°75A" in ( I3.58p . 

As a consequence on the lattice one gets (I3.4ip . and similar for the sin- 
glet axial current, for free. Of course there is a exactly conserved current 
corresponding to the symmetry (|3.58p . but it is complicated in structure and 
not ultra-local. 



The Neuberger-Dirac operator [58j is given by the expression 



D = l{l-AiA^A)-'/'} , a=Y^^, (3.60) 

A = l + s-aL>w, (3.61) 

with the Wilson Dirac operator as introduced in Section 13.491 and a pa- 
rameter < s < 1/2 that allows for some optimisation. For any A obeying 
^75 = 75^^ this operator satisfies (13.57p with a replaced by a. That is cer- 
tainly the case for More properties of the operator (13.60p are discussed 
in 60|. The best choice for the massive Dirac operator, for example, is 

Dm = {l-lam)D + m. (3.62) 



^Note that the fermion fields can be transformed independently in Euclidean space. 



Chapter 4 

The Schrodinger functional 



Nonperturbative lattice QCD is set up to deliver predictions from first princi- 
ples without uncontrolled approximations, mostly in the form of MC simula- 
tions. One ingredient on the way from numbers to physics is renormalisation. 
In order to stay on the road of first principles, also the renormalisation has 
to be done non-perturbatively. Here the Schrodinger functional (SF) has 
proved to be a powerful framework. 

The Schrodinger functional in a quantum field theory is the transition 
amplitude between field configurations at time zero and some later time. In 
an Euclidean prescription it can be represented by a functional integral 

-s 



(4.1) 

where the fields are defined on a space-time manifold with boundaries and 
obey Dirichlet boundary conditions. 

There exists a reasonable amount of literature about the Schrodinger 
functional in lattice QCD and its merits. One of the motivations to study 



the SF was, that it provides a infrared cut-off 0(1/T) to the theory [61 
where T oc L is the time extension of the lattice. Since it is a finite size 
regularisation scheme, the size of the system L is a natural scale in the theory. 
By employing finite size recursion techniques, it is thus possible to connect 
the low energy regime of the light mesons with the high energy regime of 
perturbative QCD. In this way the running of the strong coupling and the 
fundamental parameters of QCD can be computed flQiei. We can not 
give a thorough introduction to these topics here. Instead we refer the reader 
to the Les Houches lectures by Liischer [4^, which give an overview about 
the techniques. In the present chapter we discuss the renormalisability of the 
SF and on its lattice representation. 

In the continuum formulation of the Schrodinger functional the fermion 
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fields are defined in the time interval xq = [0, T]. At the boundaries they are 
subject to the conditions 

P+iIj{x) = ^{x)P- = at xo = , (4.2) 

P.ipix) = ^{x)P+ = at Xo = T. (4.3) 
The projectors occurring in these equations are defined as 

P± = i(l±7o). (4.4) 

Note that these boundary conditions are invariant under space rotations, 
parity, time reflections and charge conjugation, but not under chiral trans- 
formations (13.21) and (13.31) . 

In this chapter we concentrate on bilinear fermionic actions 

Sf = dxo j d'^xlpix) DiPix) , (4.5) 

D = 7^5m + ^ + "^0 , (4.6) 

where d = 1 or d = 3 and A{x) is the sum of bosonic fields mediating the 
interactions. In the case of QCD obviously A = 7^v4^(x), with the gauge 
field Afj_{x). In the case of the GN the four fermion interaction terms can be 
replaced by bosonic auxiliary fields in the functional integral. The fermion 
action is then also of the form (14. 5p (cf. Section ICTll . The gauge field or the 
auxiliary fields will play a spectator role most of the time. However, desirable 
properties like renormalisability and locality (in the lattice formulation) of 
the theory may depend on their details. 

In the standard formulation of the SF inhomogeneous boundary condi- 
tions are adopted and the boundary values are used as sources for the fermion 



fields at the boundary [GJJ- With the homogeneous boundary conditions of 
Eqs. (14.21 14.31) the boundary fermion fields can be defined through the non- 
zero Dirac components 

C(x) = P_^(x) , C(x) = V^(a;)P+ at Xo = , (4.7) 
C'(x) = P+^(a;) , C'(x) = V^(a;)P_ at Xo = T . (4.8) 
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This is the convention used in 

The boundary conditions of the SF may cause additional divergences. In 
quantum field theory one has learned to deal with divergences. They are 
absorbed into a "normalisation" of the fields and parameter of the theory. 
If one gets away with a redefinition of a finite number of parameters, the 
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theory is called renormalisable. It is then not obvious, that the theory with 
the boundaries stays finite. 



Symanzik argues in [63| (see also [6j| for an introduction) that the SF of 



any renormalisable quantum field theory can be rendered finite by adding a 
finite number of boundary counterterms. These are local polynomials in the 
fields and their derivatives, integrated over the boundary. They are restricted 
by the symmetries of the theory and have mass dimension d or less. 

The argument is based on explicit calculations in scalar 0^ theory, where 
one has to add two new boundary counterterms 0^ and 09o0. The expecta- 
tion could be shown to hold up to 2-loop of perturbation theory in SU(N) 



Yang-Mills theory [65|, |66|| and in QCD [45|, |67||. In the former case the 
symmetries forbid any boundary counterterms and in the later they can be 
absorbed in a multiplicative renormalisation of the quark fields at the bound- 
ary . Although a proof is still missing, there is little doubt that the SF of 
QCD is renormalisable, given the success of the method in non-perturbative 
computations [68|. 

Discretising the Schrodinger functional of a given quantum field theory 
involves then two questions. Boundary conditions like Eqs. (14.21 14.31) make 
only sense when imposed on smooth functions, but the lattice is discrete by 
definition and the lattice fields are only defined at the sites of the lattice. So 
the question arises how the continuum boundary conditions are represented 
on the lattice. Furthermore the lattice breaks some continuum symmetries, 
like continuous rotations, translations and chiral transformations. This may 
give rise not only to new terms in the bulk of the lattice, but also to additional 
terms at the boundaries. 

In Section 14.11 we address the first question with an heuristic argument 
to make plausible, that the Schrodinger functional boundary conditions arise 
naturally in the continuum limit and need no special adjustments of the 
lattice action. Then in Section 14.21 we get explicit and present formulae for 
the well studied lattice SF with Wilson fermions. We present a rather recent 



proposal [32| for an overlap Dirac operator in the SF in Section [43l There is 



another proposal for a SF with Ginsparg- Wilson fermions by Taniguchi [69 
There the boundary conditions are realised through an orbifold projection. 
However, there are some technical difficulties, i.e the fermion determinant 
has a phase and fermions masses can not be introduced straightforwardly. 



Therefore we consider only the proposal of Ref. [33|. At the end of this 



Chapter we come back to the question of boundary counter terms. 
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4.1 Lattices with boundaries 

Discretising a field theory involves some freedom in the definition of the 
lattice action. For example the differential operator can be replaced by the 
forward finite difference ad^ip{x) = ip{x + aiJi)—ip{x) or by the averaged finite 
difference ad^ip^x) = ip^x + a/x) — ipi^x — a/i). That both choices lead to the 
same continuum theory is assured by the local nature of the resulting actions. 
In the continuum limit the relevant length scale of the theory becomes much 
larger than the lattice spacing a and the microscopic (length scale a) details 
become irrelevant. One says the two actions belong to the same universality 
class. 

These classes are characterised by global properties like dimensionality 
and symmetries. In statistical mechanics universal behaviour occurs in the 
vicinity of the critical lines and there the concept of universality has been 



extended to systems with boundaries (see [70!| for a review). This means for 
the discretisation of field theories that boundary conditions imposed in the 
continuum theory, together with the dimensionality and the symmetries of 
the theory, define an universality class. In general requiring locality, symme- 
tries and power-counting reduce the number of possible boundary conditions 
and therefore the number of universality classes. 

The line of argument followed here is borrowed from Section 3 in Ref. 
(3^ and we also start with the simple example of a free scalar field. 



4.1.1 Free scalar field 

We consider a free scalar field in the half-space Xq < 0. We do not specify 
the boundary condition at Xq = 0, but rather explore where we are lead to 
in the continuum limit starting from different lattice actions. The lattice 
fields (/)(x) are defined at the sites of a hypercubic lattice with spacing a. 
Only the fields at xq = a, 2a, 3a, . . . are dynamical degrees of freedom and 
are integrated over in the functional integral. A possible lattice action reads 

S = a^J2Yl H9M^)9,<Pi^) + m'<l>{xy} , (4.9) 

where 9^ is the forward difference operator in direction. Note that the 
action depends only on the dynamical degrees of freedom. 

In order to calculate the propagator we write Eq. (14.91) as a quadratic 
form in (j) 

S = a''J2Yl + + PMx) (4-10) 

xo>a X 
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with the backward difference operator 9^. Replacing the symmetric nearest 
neighbour interaction of Eq. (14.90 by the one in Eq. (14.101) produces some 
mismatched terms near the boundary. This is cured by the inclusion of the 
boundary term 

P = -^4o,a5*. (4.11) 

The defining equation of the propagator G(x, y) follows directly from the 
action (liHHl 

i-d;d^ + m^ + P)Gix,y)=a-^5,,y, xo,yo>a. (4.12) 

Due to translation invariance in the space directions the propagator can be 
calculated in a time-momentum representation 



G{x, y) = e'^^-~y^ G{xo, yo, p) • (4.13) 



T/a 

This is done by first determining the eigenfunctions of —(9^(9^ + that are 
annihilated by P 

e^P" cos ((a;o - f )po) ■ (4.14) 
It is then easy to write down the propagator 

/2vr/a 
~^ m? +p'^ ^ ^^^^^ ^ ^^^^^ ■ ^^'^^^ 

In the continuum limit the integral can be done and the result is 

G'cont(x, y) = j (e-(p)l-o-.ol + e-.(p)(.o+.o)) ^ (4.16) 



e(p) = Vm2 + p2. (4.17) 

Examining the behaviour at the boundary one infers that the propagator 
satisfies Neumann boundary conditions in the continuum limit 

5oG'cont(a;,l/)L^^o = 0. (4.18) 

Now we slightly modify the action by adding a boundary term 

S^S + a'Y,^^i^f _ ■ (4-19) 



X 



Note that the powers of a are such that c > is dimensionless. The propa- 
gator is still given by (I4.12p but now P has an additional term 
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The eigenfunctions that are annihilated by P are a bit more complicated 

e^P'' sin ((xo - a)pQ + ^{po)) , (4.21) 

^(po) = arctan { f\^^'^ ) . (4.22) 

Vcos (po) +C-IJ 

The propagator is similar to (I4.15p but with the eigenfunctions f l4.2ip . Again 
the integral can be performed in the continuum limit and the result is 

Wf (e-^(P)'--^"' - e-^(P)(-+^°)) . (4.23) 

This propagator satisfies Dirichlet boundary conditions 

G':ont(^,Z/)L.„=o = 0- (4-24) 

A small change in the action led to a totally different class of boundary 
conditions. This can be made a bit more transparent by looking at the field 
equation 

(5 S 

{vix)(Piy)) = a-%,y , ri{x) = . (4.25) 

From the action in the form (14.191) we derive 

r]{x) = {-d;d^ + m^} (P{x) , xo>a. (4.26) 

Thus in the bulk of the lattice we find the Klein-Gordon equation. But at 
the boundary Xq = a we find 



r]{x) = 4 --do 0(x) + {-dldk + rr?} 0(x) . (4.27) 



The negative powers of a in front of the first two terms are due to the fact 
that they are supported only at the boundary. Taking the limit a — >■ the 
first term dominates for all c > 0. Thus in the continuum limit the field 
equation at the boundary implies Dirichlet boundary conditions at Xq = 0. 
In the somehow special case c = the second term dominates and the field 
equation implies Neumann boundary conditions. 



4.1.2 The SF universality classes 

In the example of the last section no boundary conditions were imposed 
onto the lattice fields. The boundary conditions were encoded in the lattice 
action and emerged when we took the limit a — > 0. But we found two distinct 
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classes of boundary conditions. One that is generic, in the sense that it is 
found for a wide range of actions (all c > 0). And one that is sensitive to 
small perturbations of the action. 

Clearly a free scalar field is a trivial example. In the case of interacting 
theories or more complicated actions the analysis of the field equations will 
not be so transparent. But boundary conditions that respect locality will 
always be of the form 

^(^)Lo=o = 0' (4-28) 

where 0{x) is a linear combination of local fields and their derivatives with 
the appropriate symmetry properties. 

Inspired by the above example it is plausible that the generic boundary 
conditions, that are stable under perturbations of the lattice action, are those 
imposed on the fields with the lowest dimension. All other possible boundary 
conditions will require some tuning of the lattice action, unless there are 
symmetries that protect them. 

We want to apply this argument now to the case of fermions in two (GN) 
and four dimensions (QCD). In both cases the fermions are represented on the 
lattice by Dirac spinors ip{x) at each lattice site. Both theories are asymptot- 
ically free, assuring that the scaling dimension of local fields is equal to their 
engineering dimension (cf. Section [2TT1) . Thus the fields of lowest dimension 
are the fermion fields themselves and the generic boundary conditions are of 
the from 

B^{x)l^^^ = 0, (4.29) 

where -B is a constant matrix with Dirac, fiavour and, in the case of QCD, 
colour indices. Boundary conditions at a later time xq = T are then linked 
to those at Xo = by time refiection symmetry and the conditions for the 
anti- fermion field are given by charge conjugation symmetry. 

The possible boundary conditions are further restricted by the lattice 
symmetries. If the lattice theory is invariant under gauge and fiavour trans- 
formations, cubic rotations and parity, so should be the boundary conditions. 
And finally the matrix B cannot have full rank since the Dirac equation is a 
first order differential equation and the two conditions at Xq = and Xq = T 
would imply a vanishing fermion propagator. Up to a constant factor there 
are only two matrices satisfying all these properties 

B = P_ and B = P+ . (4.30) 

The homogeneous Dirichlet boundary conditions of the Schrodinger func- 
tional are therefore the generic ones and do not need any fine tuning or par- 
ticular adjustment of the lattice action. There are two universality classes of 
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lattice theories which differ by the sign in the boundary condition 

P±^(x)L^^o = 0. (4.31) 

Since the two are connected by a finite chiral transformation, the difference 
matters only for non-vanishing fermion masses. In this case the sign can be 
determined by inspection of the free propagator. 



4.2 Free Wilson fermions 

In the continuum the Schrodinger functional can be represented by a func- 
tional integral 

Z = j DipD^e-^^^^^'^^ , (4.32) 

where the fields obey Dirichlet boundary conditions (l4.2H4.3p and the action 
may be given by (14. 5p . For simplicity we omit any interaction (gauge or four 
fermion) and introduce a d + 1 dimensional lattice with spacing a and label 
the sites by integer multiples x^/a G Z , = 0, . . . ,d. The space directions 
are taken to be periodic with length L 

V^(x + Lk) = ij{x) , (4.33) 

where k is an unit length vector in direction k = 1, . . . ,d. The fermion fields 
at times xq = a, 2a, . . . ,T — a are the dynamical degrees of freedom (the 
fields that are integrated over in the functional integral). It is convenient to 
assume that the fermion fields are defined at all other values of xq as well, 
but that they are zero there. The free lattice action takes then the familiar 
form 

So = a''J2^{x)Dmij{x), (4.34) 

X 

where Dm is some discretisation of the massive Dirac operator. 

4.2.1 Dirac operator and propagator 

In the present section we consider 

Dm = Dv^ + mo, (4.35) 
where the Wilson Dirac operator is defined as (cf. Section [3^ 

= + d,) - ad;d,} , (4.36) 
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in the range < xq < T. At all other times the target field x = Dmip is set 
to zerc0 



xo<0 



= Dm^{x) 



xo>T 



(4.37) 



Note that we set r = 1 and that we introduce factors e*"^*^/^ in the spacial 
lattice difference operators (which is equivalent with periodic boundary con- 
ditions with a phase, see Section [2T2ll . For 9 = this is the Dirac operator 
introduced in 

Thus the Dirac operator can be considered as a linear mapping in the 
space of fermion fields that vanish at the boundaries. The propagator S{x, y) 
is defined through 



with boundary values 



1 



a 



d+l "^'2/ ' 



0<x,y <T , 



P+S{x,y) 



xo=0 



P-S{x,y) 



0. 



xo=T 



(4.38) 



(4.39) 



Since the operator 75-Dm is hermitian the propagator has the property 

S{x,yy = ^,S{y,x)^,. (4.40) 

In the free theory it is possible to derive an explicit expression for the 
propagator in a time-momentum representation |7l|. An elegant form is 



S{x,y) = {D\^ + mo)G{x,y) , 0<xo,yo<T, 
where G{x, y) is defined through 
G{x,y) =L-'^ 5^{-2zp+A(p+)i?(p+)}-' e'P(--y) 



(4.41) 



X I (M(p+) - ip+) e-^(p)l^«-2^ol + (M(p+) + ip+) e-^(p)(22^-No-i'ol) 
-(M(p+) + z7opJ) e-"(P)("o+^") - (M(p+) - nopj) e-^ip)i^T-x,-y,) 

(4.42) 

All undefined functions and notations in this expression are introduced in 
Appendix IA.2[ In particular the momenta are given by eq. IA.19I and 



^The Dirac operator maps the space of fermion fields that are defined at all xq, but are 
zero aX xq < a and xq > T — a, into itself. 
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4.3 Ginsparg- Wilson fermions 

The Wilson lattice Dirac operator violates chiral symmetry explicitly and 
leads to computational difficulties such as additive mass and multiplicative 
current renormalisation. Therefore it is desirable to have a lattice operator 
with better chiral properties. In Section 13.31 we learned that it is not pos- 
sible to formulate a lattice theory of fermions with a continuum like chiral 
symmetry. The way out in Section [331 was to ease the restriction to an exact 
continuum like chiral symmetry in favour of a lattice chiral symmetry and a 
lattice Dirac operator given as the solution to the Ginsparg- Wilson relation 

-f5D + D-f5 = aD-f5D. (4.43) 

In the limit a ^ this relation becomes the known anticommutation relation 
and the lattice chiral transformations become the continuum ones. 

In the presence of the boundaries the same approach immediately leads 
to inconsistencies. As stated in the introduction in the continuum the SF 
boundary conditions break chiral symmetry. This can be seen considering 
the solution to the massless Dirac equation 

Dip{x)=0, (4.44) 

which has to obey the boundary conditions (l4.2H4.3p . Note that the mass- 
less Dirac operator anti-commutes with 75. Therefore, given the propagator 
D S{x, y) = 6{x — y) the sum 75 S{x, y) + S{x, y) 75 is a solution of the Dirac 
equation for all y. And since the propagator obeys the boundary conditions 
f l4.2H4.3l) this solution can also be obtain from the boundary values at xq = 
and xq =T (actually the non-zero components there) 

l^SiyX.y) + S'(x,y)75 = 

/ d''zS{x,z)^,S{z,y)+ [ d''zS{x,z)j,S{z,y). (4.45) 

Jxo=0 Jxo=T 

Thus the continuum propagator anti-commutes with 75 up to boundary 
terms. In QCD with more than one massless quark this leads to a non- 
singlet chiral Ward identity with a unit mass term at the boundaries (s^ . 

Now consider a lattice Dirac operator satisfying (14.431) in the presence of 
the boundaries. Then the lattice propagator anti-commutes with 75 

75 5'(x, y) + 5 {x, y) 75 = ^ 6^^y-f5 , (4.46) 

for any finite separation x — y and any finite a. Thus this theory can not have 
the right continuum limit, i.e. the one where the propagator obeys fl4.45p . 



4.3 Ginsparg- Wilson fermions 
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4.3.1 Modified Neuberger-Dirac operator 

As explained in Section [4.1.2l the SF boundary conditions form a universality 
class and need no special adjustment of the lattice action. This means in 
particular that there may be many possible lattice operators leading to the 



right continuum limit. In Ref. |32| it is proposed to allow for an additional 



term on the right hand side of Eq. (I4.43P 

+ = a^75^ + Ab, (4.47) 

where is supported in the vicinity of the boundaries and decays expo- 
nentionally with the distance to them. 

The operator introduced in the same reference 

D^ = ^{1-'^{U + U-)} , (4.48) 

U = A {A^A + caP) , U- = 75f/^75 , a = , (4.49) 

is a modification of the Neuberger-Dirac operator (13.601) . Here A is essentially 
the Wilson-Dirac operator in the presence of the boundaries f l4.361l437ll 

A=l + s- aD^ . (4.50) 

The parameters s and c can be used to optimise numerical computations. 
The modification is due to the boundary operator 

P^(X) = ^ [5.,,aP-i^{x)l^^^ + 5,„T-aP+i^{x)l^^^_^] • (4.51) 

The construction is such that inherits the transformation properties 
of the Wilson-Dirac operator under cubic rotations, parity, time-reflections 
and charge conjugation. And having the combination U + f/~ renders 75-Dn 
hermitian. Although we will be mainly concerned with massless fermions the 
extension to massive fermions is simple ^] and given by 

= {I - \amQ)D^ + mo . (4.52) 

In Section [4.3.31 we will show that the operator fl4.48p indeed obeys fl4.47p 
with a replaced by a. 
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4.3.2 Free theory 

For our perturbative computation we need the Dirac operator and the fermion 
propagator in the free theory (no gauge field, no four fermion interactions). 
To obtain an exphcit expression for the free propagator like in the Wilson 
case might be possible but is much more difficult. Instead it is computed 
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numerically using well established methods 

Nevertheless, the operator under the square root in (14.491) can be worked 
out in the time-momentum representation and it is reassuring that for the 
interesting range of the parameters s and c its eigenvalues are positive. This 
can be seen as follows. In the free theory the operator under the square root 
explicitly reads 

A^A + caP = (1 + + so" d;d^ + y Yl d;d^dtd, + (c - l)aP , (4.53) 

and acts on the fermion fields in the presence of the boundaries (cf. Section 
I4.2.1I) . It is hermitian and therefore has real eigenvalues. Due to translation 
invariance in space the spatial eigenfunctions are plane waves 

«^(^) = i;?E ^p(^o)' (4.54) 
p 

with spatial momenta pk in the range (IA.33I1 . In the following we restrict 
ourself to a definite momentum p. Then the Dirac components of Wp{xo) are 
ordinary functions of Xq. Therefore we drop the subscript in the eigenvalue 
equation 

(— g a^dodo + m + (c — l)aP) w{xo) = Xw{xo) , (4.55) 
where q and m are short hand for 

k<l 

Because of the projectors in P, eqs. (14.551) are two coupled equations for 
the plus and minus components in w{xo) = P^w^{xo) + P_ w_(xo). In the 
special case c = 1 they have to be identical and one easily finds the solutions 

TITT 

sm{poXo), pq = —, n = 1,2, ... ,T/a - 1 , (4.57) 
and the correspondent eigenvalues 

\ = {l + sY- saY + ^YpIp- = QPl + m. (4.58) 



4.3 Ginsparg- Wilson fermions 
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For arbitrary c the solutions to (I4.55p are 

P_ sin {poxo + 6) + P+ sin {po{T - xq) + h) , (4.59) 

h = - arctan ( sm(QPo) — ^ \ 
\^^ + cos(apo) y 

and the allowed values of po are given by the solutions of the equation 

tan poT = r . 4.61 

^ + cos(apo) 

The solutions of this equation are all real for non-negative q/{c— 1). Since 
the eigenfunctions are odd functions of po it is sufficient to stick to the T/a — 
1 solutions in the interval [0,7r/a). In any case the eigenvalues f l4.58p are 
bounded from below by (1 — |s|)^. 

For g/(c — 1) < one finds a pure imaginary solution pg = ik, k > 0. 
This solution enters the eigenvalue (14.581) through the product qpl and might 
lead to zero modes and negative eigenvalues. However, for c > 1 imaginary 
solutions exist only for g < and thus the product q (Po)^ is positive and the 
eigenvalues (I4.58P are as well bounded from below by (1 — |s|)^. 

Since the operator under the square root is bounded from below we can 



adopt the argument in [73|, using expansion in Legendre polynomials, to 
conclude that in the free theory the locality of the Dirac operator is 
guaranteed for all |s| < 1 and c > 1. The eigenfunctions (14.591) can be 
orthonormalised and used to write down an analytical expression for the 
kernel -Dn(x, y) of the Dirac operator. But the evaluation of Df^{x, y) in this 
way would be very expensive, since it involves a sum over momenta po which 
in turn are determined for each set of parameter values by the roots of (14.611) . 

Even more desirable would be an explicit expression for the propagator. 
With the eigenfunctions of the operator under the square root at hand one 
would hope to find the eigenfunctions of the hermitian operator Dn-DJj and 
write the propagator as 



S{x,y) 



D 



{x,y). (4.62) 



But the eigenvalue problem of -Dn-DJj is much more complicated than the 
one treated above and we do not attempt to solve it here. Instead the prop- 
agator is computed numerically using established methods. This includes a 
polynomial approximation of the inverse square root in the Dirac operator 
and solving the linear equation 

D^iP = V, (4.63) 
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Figure 4.1: The deviation from the Ginsparg- Wilson relation a As in Eq. I|4.64p is 
localised at the boundaries with tails that decrease exponentially with the distance from 
the boundaries. The plot is for a 16 x 32 lattice, 9 = 0.5 and xi — yi — L/2. 



for xjj with the approximate Dirac operator and appropriate sources r?. The 



precision can be controlled throughout all steps of the computation [72|. In 
this way the whole propagator can be evaluated with a given accuracy. In 
practice we also Fourier transform to the space component to obtain the 
time-momentum representation (see Section r5.5.2l) . 

4.3.3 Chiral properties 

As a check, it was tested numerically whether the operator (14.48p is a solution 
to 

75D + D75 = aD75D + As, (4.64) 

with 

||aAB(x,y)|| ^e-'-'/' + e-'C^-")/", /t > , (4.65) 

for large distances r and T — r from the boundaries. For simplicity we define 
r as the distance from the boundary at Xq = along the time direction as 

T = — ^ — . (4.66) 

In Fig. [m] aAsix.y) is plotted for a 16 x 32 lattice and = 0.5. The 
deviations from the Ginsparg- Wilson relation decay exponentially and from 
the "effective mass" plot Fig. 14.21 the rate seems to approach k = \/2 for 
T ^ 00. It should be possible to compute AB{x,y) analytically to some 
extent, but we did not attempt to. 



4.4 The generating functional 
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Figure 4.2: "Effective mass" plot for aAsix^y). For large distances from the boundary 
aAB{x,y) is compatible with l|4.65p with k w \/2 (dotted Hne). The plot is for the same 
parameters as Fig. 14.11 



4.4 The generating functional 

We want to compute expectation values of polynomials O in the fermion and 
anti-fermion bulk and boundary fields. A possible lattice representation of 
the boundary fields isl 

C(x) = P_^(x)| , C(x)=^(x) P+, (4.67) 



I xo=a 



xo=a 



C'(x)=P+^(a;) , C(x)=^(x) 

xo=T—a 



P_ . (4.68) 

xo=T-a 



With this choice it is enough to introduce sources for the fermion and anti- 
fermion fields r]{x) and r/(x) in the interior of the lattice < xq < T. The 
generating functional is then 

Z\rj,ri] = J DipDip exp <|— 5*0 [■?/', ^/^j + a"^^^ [ilj{x)r]{x) +fj{x)ilj{x)] 

(4.69) 

As said before only the fields at < Xq < T are integrated over in this 
functional integral. We can perform this integration and obtain 

\nZ[r], T]] = a^('^+^) ^ fj(x) S{x, y) r]{y) + const. . (4.70) 

x,y 



^In QCD one has to include link variables into this definition of the boundary fields. 
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The Schrodinger functional 



Thus the generating functional is an exponential of a quadratic expression in 
the sources. Replacing the fields in O by functional derivatives 

we may write the expectation value as 

{0) = {]-0z\ . (4.72) 

L ^ J rj,r?=0 

They are given as the sum of all Wick contractions. Below we list the basic 
contractions. 



[%lj{x)^{y)] 


= S{x,y) , 






(4.73) 


[V^(^)C(y)] ^ 


= S{x]a,yi)P+, 






(4.74) 




= S{x;T - a,yi) 






(4.75) 


[C(x)^(z/)] = 


= P-S{a,xi;y) , 






(4.76) 


[C'(x)^(y)] 


= P+S'(T — a, xi, 


~y) , 




(4.77) 


[C(x)c'(y)] 


= P^S{a, xi] T — 




)P-, 


(4.78) 


[C'(x)C(y)] 


= P+S(T — a, xi] 




)P+, 


(4.79) 


[C(x)C(y)] = 


= P-S{a,xi; a,yi 


)P+, 




(4.80) 


[C'(x)C'(y)] 


= P+S(T — a, xi 


;T- 


a,yi)P_ . 


(4.81) 



4.5 Boundary counter terms 

The boundary conditions of the SF may give rise to new counter terms defined 
on the boundary, in the sense that the associated bare coefficients are needed 
to absorb infinities on the way to the continuum limit. Relevant operators 
or composite fields living at the boundary are those which have dimension d 
or less. Thus objects like ipTip with T = {1,70,71,75} can appear. 

The identity and 70 can be written as P+ + P_ and P+ — P respec- 
tively. Hence the corresponding terms are proportional to the boundary 
fields (14 .714. 81) in the continuum or (I4.67|4.68p on the lattice. The terms 
with 7i and 75 violate parity and are therefore not present. 

Thus it is enough in the course of renormalisation to introduce a renor- 
malisation factor for all boundary fields 

C(x)^ZcC(x), C(x)^^cC(x), (4.82) 

and equivalently for ( . 



Chapter 5 



Self-coupled fermions in two 
dimensions 



5.1 Four fermion operators 

In two dimensions fermion fields have mass dimension 1/2 and a local four 
fermion operator 

'^ai'^Pj^lklpSl , (5.1) 

has a dimensionless coupling in the action. From the point of view of dimen- 
sional analysis such a interaction term is renormalisable. This means that if 
such a four fermion interaction is added to a theory that is renormalisable, it 
stays so. All new divergences can be absorbed into a redefinition of the four 
fermion coupling. 

In the literature one finds several two-dimensional fermion models with 



different symmetry and interaction content [19|, l20|, l23||. This is because 
there are a lot of possible ways to contract the indices of four fermion fields. 
But in two dimensions there are also several relations between the possible 
contractions. In the following we discuss this in some detail. 

The most general operator is an arbitrary contraction of the Dirac (Greek 
letters a, (3, . . . ) and fiavour (Latin letters ) indices in eq. (|5.1I) . Since 

we want a theory that is as similar to QCD as it can be in two dimensions, it 
certainly should be Lorentz invariant, even under parity and, in the massless 
case, have an U{N) fiavour symmetry. 

Lorentz invariance strongly constrains the possible contractions of Dirac 
indices. Since eq. (15.11) must be a Lorentz scalar it must be a product of two 
scalars, two pseudo-scalars, two vectors or two axial-vectors. Consider the 



^This is meant in the continuum. At the end of the day we are interested in continuum 
QCD. 
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case of two scalar^ 

i^ai^aj^pki^pi ■ (5.2) 
Invariance under U{N) transformations 

where U^U = UU^ = I , (5.4) 
allows for the following two flavour contractions 

^m^m^/3j^/3j ) ^ ai^aj^ pji' (ii ■ (5.5) 

The second contraction can be expanded in a basis of N x N matrices (Ap- 
pendix [AAl) 

a 

where we suppressed all subscripts on the left hand side. The first term is 
proportional to the first one in (15. 5p and the sum is over the generators A'* 
of SU{N), which act on the fiavour indices of ip and ip. 

Therefore the most general four fermion operator consistent with the 
symmetries can be expanded in a basis of eight different contractions. One 
half of them are products of flavour singlet bilinear operators 



Oss = {m\ 

Opp = (V'75V;)^ 

Ovv = E^(V'7mV') 



2 



(5.7) 



Oaa = E;.(^7m75^)^ , 
and the other half are products of flavour vector bilinear opeartors 

Oi^y = E,,.(^7.A'^^)^ ^'-'^ 

But in two dimensions these operators are not independent. Suppose 
there is only a single fermion (A^ = !)■ From the above list one would expect 
four different terms (no flavour vector operators for a single fermion). But 
there are only four independent field components at each space-time point. 
Thus there is only one local four fermion operator. 



^Repeated indices are summed over if not indicated otherwise. 



5.2 Chiral symmetry 
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For > 1 the number of independent field components is no restriction. 
Nevertheless there are relations among the operators above that reduce the 
number of really independent ones to three. Because of a peculiarity of the 
7-matrices in two dimensions (7^75 = ie^vlv) there is no difference between 
vector and axial- vector and thus 

OaA = -OVV , O'aa = -O'yy . (5.9) 

More dependencies are due to Fierz identities. Fierz identities connect 
products of Dirac bilinear forms by rearranging the order of the Dirac spinors 
(see Appendix IA.3.11 ). For the two flavour contractions of (15. 5p . but with 
general Dirac structure, we find 

(V^,r^,) (v^/^,) = -\Y. tr(r, r r, r)(v^,r,v^,) (v^/,v,) , (s.io) 

where F, F/ G {F5 = 1 , Fp = 75 , Fy = 7o , 7i}- Using also (15. 6p this yields 
three identities relating the fiavour-singlet and the flavour-vector operators 

O'ss = -{l + 2/N)0ss-0pp-0vv (5.11) 
O'pp = -0ss-{l + 2/N)0pp + 0vv (5.12) 
O'yy = -20ss + 20pp-2/NOyy. (5.13) 

The number of independent operators has been reduced from eight to 
three. One possible choice would be 

Oss.Opp.Ovv (5.14) 

But any other combination of three operators from (|5.7I) and (j5.8l) is equally 
good. 

5.2 Chiral symmetry 

The eight operators of ( I5.7p and (15.81) . of which only a set of three is indepen- 
dent, are invariant under U{N) transformations (15.31) by construction. Since 
massless QCD at the classical level has a global chiral U{N) symmetry we 
investigate here the transformation properties of these operators under the 
infinitesimal transformations (13.241) and (13.261) . In Section [34] these were in- 
troduced as local transformations. But since they will act on operators that 
are local and contain no derivatives of fields the outcome is also applicable 
for the global case. 
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For the above mentioned symmetry transformations the change 60 of the 
operator O linear in the infinitesimal parameter u is 

O ^O + 6O + 0{uj^). (5.15) 

If SO vanishes the operator is invariant under finite transformations that can 
be build up from the infinitesimal ones. 



Axial f/(l) transformation Consider first the axial U{1) transformation. 
Because of the 7-matrices in Ovv and Oyy they are invariant like the kinetic 
term. 

6Ovv = 0, 6O'yy = 0. (5.16) 

But SOss and 60pp do not vanish. Using the projectors Pr^l = |(1 ± 75) 
we write 

Oss = {Pr + Pl) i^Y , Opp = (^ {Pr - Pl) V^)' . (5.17) 
Now it is easy to see that 

60 SS = 2tUJA{^{Pn-PL)^){^{PR + PL)^), 
60pp = 2lU0A{^{PR + PL)i'){i'{PR-PL)ij). 

Therefore the difference of the two operators is invariant 

6{0ss-0pp) = Q, (5.18) 
and similar for the primed operators 

5{0'ss - O'pp) = Q . (5.19) 

Axial SU{N) transformation The axial SU{N) transformation (13.261) 
affects Dirac and fiavour indices. However, Oyy has a trivial fiavour structure 
and it is easy to see that it is invariant 

6Oyy = 0. (5.20) 

Computing 60yy the commutator [A", A^] = 2if°-^'^ with the structure con- 
stants /"'"^ appear (cf Appendix lA. 3. 2^ 

60'yy = 2zw^(V^75A'^7^A^V> + V^7^_A''A'^75V>)(V^7^A''7/>) 

= -2zu;^(V^_7^75[A", A"] ^) (^7^^' V^) 

= Mr^^(^7M75A'=^)(^7MA'^) 

= 4^a;y'^'"={(V'7iA^^) (V'ToA^'V') - (V'ToA^V') (V'TiA^V')} ■ 



5.3 Lattice chiral Gross-Neveu model 
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In the last step we used the definition of 75 and the Cliff'ord algebra (cf. 
Appendix lA.l.ll) . Since the structure constants are totally anti-symmetric, 
the whole expression is symmetric in the indices b and c and hence does not 
vanish for general a;^ 

50'yy^O. (5.21) 

The result for Oyv and O'yy can be used together with fl5.11f[5T3l l to infer 

6{Oss-Opp)y^0, and 6{0'ss - O'pp) . (5.22) 



5.3 Lattice chiral Gross-Neveu model 

The continuum chiral Gross-Neveu model in Euclidean space-time is given 
by the action 

-^CGN = / d'x {i^i,d,tp- \g\Oss - Opp) - \glOvv} ■ (5.23) 

From the analysis of the last two sections we know that this is the most 
general action with chiral U{1) x U{1) and SU{N) fiavour symmetry. But 
we also learned that this form is not the only possibility. Using the identity 
(15.131) the action (15.23p may equally be written as 

ScGN = J d'x {i^i,d,^ + Ig'O'yy - i^Oyy} , (5.24) 

with 5l = gl- g'^/N . (5.25) 

The coupling 5y is believed to have an exactly vanishing beta-function. This 
is connected to the decoupling of the U{1) part of correlation functions that is 
derived by formal manipulations of the continuum path integral |74l.l75|. The 
vanishing of the beta-function has been proved up to two-loop perturbation 
theory (S). 

Passing now over to the lattice inevitably breaks part of the continuum 
symmetries. Generally this leads to more possibilities for the mixing of oper- 
ators under renormalisation and to additonal parameters in the action. For 
example, an operator that is multiplicative renormalisable in the continuum 
may lose this property on the latticel^ Being lattice artefact^! one has to 
assure that these effects disappear in the continuum limit. As we will see. 



^This also occurs in other regularisation schemes. In the two- loop computation of [22| 
using dimensional regularization so called "evanescent operators" appear at an intermediate 
state of the calculation. 

^This is not true if there is an anomaly that breaks the symmetry. 
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this may happen automatically, but in some cases additional parameters in 
the action have to be tuned. 

At first the theory is set up on a hypercubic lattice with periodic boundary 
conditions. The introduction of boundaries as in the Schrodinger functional 
involves additional problems that are addressed in a separate section (Section 
EE). 

In the action we have to include terms that are forbidden by the contin- 
uum symmetries, but are allowed by the less restricting lattice symmetries. 
By dimensional analysis such operators with mass dimension n will have a 
factor a"~^ in front (where D = d + 1 is the dimension of space-time and 
should not be confused with the Dirac operator). Operators with n > D can 
safely be neglected since they automatically disappear for a 00 But the 
operators with mass dimension n < D are needed and lead to an additional 
renormalisation. 

Obviously (Euclidean) Lorentz symmetry 0(2) is broken on the lattice. 
The symmetry group of the hypercubic lattice, the hypercubic group, is a 
subgroup of 0(2) and contains rotations by 7r/2 and reflections. Luckily, the 
operators allowed by the hypercubic group are the ones already present in 
(I5.23P or (I5.24P and operators with mass dimension n > D. So, ignoring all 
other symmetries for a moment, the efl'ects of the broken Lorentz symmetry 
in the action automatically disappear in the continuum limit and we may use 

^CGN = {i,Di,- \g\Oss - Opp) - \glOvv] , (5.26) 

X 

with the lattice Dirac operator D as the lattice action. The fermion fields 
ip{x) and ip{x) are defined at the sites of a hypercubic lattice with periodic 
boundary conditions and lattice spacing a (c.f Section [2?2ll . Note that there 
may be mixing amongst operators that cannot be neglected and that thus 
has to be taken into account when calculating expectation values of such 
operators. 

Depending on the exact definition of D there are further broken or mod- 
ified symmetries that might lead to a modification of f l5.26p . The Wilson- 
Dirac operator explicitly breaks chiral symmetry and leads to an additive 
mass renormalisation. For staggered fermions the broken fiavour symmetry 
causes new mixings. The lattice chiral symmetry associated with Ginsparg- 
Wilson fermions forbids additional operators of mass dimension n < D in 
the action. In the next two sections we discuss Wilson and Ginsparg- Wilson 
fermions. 



^The n — D + 1 operators may be considered if one wants to achieve 0(a) improvement. 



5.3 Lattice chiral Gross-Neveu model 
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5.3.1 Wilson fermions 

As we have seen in Section 13.41 the Wilson-Dirac operator explicitly breaks 
chiral symmetry and a fine tuning of the bare mass is needed to restore chiral 
symmetry. The action should therefore contain all terms allowed by the 
remaining U{N) flavour symmetry. In Section ISTTl a basis for U{N) invariant 
four fermion operators was established. Due to dependencies among the 
diverse operators it is enough to choose as set of three operators out of f l5.7p 
and (15. Sp . In analogy to (|5.26l) one possibility is 

<S'cGN,W = «^ 1^ (^W + "^o) i' 

X 

- jg\Oss - Opp) - i^Opp - Ig'yOvv} , (5.27) 

where we added a chiral symmetry breaking mass term and the coupling of 
Oss and Opp are no longer related. 

There are many different choices one could make. For example, one can 
invert eqs. ( 15.1ip -( l5T3ll to obtain 

5'cGN,w = X] + "^o) 

X 

- \g'\0'ss - O'pp) - \5'lO'pp - ylO'yy] , (5.28) 
where the primed couplings are related to the unprimed ones 

■[{2/N)g^-5l-2gl], (5.29) 
^P = i(iV^{2(ViV-l)^p}' (5-30) 

= 4(iV^ - 1) { - + + ■ (5.31) 

The traceless A-matrices in the four fermion operators in (15.281) can help to 
reduce the number of non-zero diagrams in perturbative calculations. (This 
will be utilized in Section [6]) Yet another version is obtained by decomposing 
the interaction as in (15.240 

'S'cGN,w = X] (^w + mQ)i) + jg^Oyy - \5pOpp - \5lOyy] , 

(5.32) 

with 5l = gl- g'^/N . (5.33) 
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As pointed out in Section [3^ the bare mass mo lias to be tuned to a non- 
zero value in order to restore chiral symmetry. In the same way 6p has to be 
tuned since the value 6p = 0is not distinguished by a greater symmetry. A 
whole section (Section [6l) is devoted to the restoration of chiral symmetry. 

5.3.2 Ginsparg- Wilson fermions 

A Dirac operator satisfying the Ginsparg- Wilson relation f l3.57p implicates 
the lattice chiral symmetry (13.581) . This symmetry becomes the familiar 
chiral symmetry in the continuum limit and forbids a mass term. However, 
at finite lattice spacing the invariant four fermion operators are diff'erent from 
the continuum ones. 

The transformations (I3.58p may be written as 

tp ^ip + e%ip , ip ^ip + eip'j^ , (5.34) 

with 75 = 75(1 — aD) and 75 = 1 . (5.35) 

This looks like the usual chiral transformations except for the %. Using 
f l3.58p one proofs the identity 

(l-fD)75 = 75(l-f/^). (5.36) 
Which means that (1 — ^D)ip transforms under fl5.34p like ip under fl3.58p 

(1 _ |Z))^ in)ij + 675(1 - ID)ij . (5.37) 

Therefore the chirally invariant operators of Section 15.21 are invariant under 
the lattice symmetry after replacing ip ^ {1 — f-D)^'- In particular, if we 
define _ 

Opp = (^75(1 - ^D)^f , (5.38) 

then the operators invariant under 15.21 are 

6dvv = 0, 6{dss-0pp) = 0- (5.39) 

The treatment and result in the case of the operators which are products of 
flavour vector bilinear operators (15.81) is equivalent. Note that the operators 
f l5.38p are not ultra-local any more. Nevertheless the identities (15.91) and 
(I5.1ip - (l5.13p hold also for these operators. But in the derivation in Section 
15.11 one has to replaces ip ^ ip = {1 — ^D)ip. 



5.4 The discrete Gross-Neveu model 
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Now we gathered all prerequisites to write down the lattice action for 
Ginsparg- Wilson fermions in the form 

^CGN,GW = a' 5^ {i^D^lj- yidss - Opp) - \gldvv] , (5.40) 



or 



5cGN,GW = a'$^ [^D^ + \g^d'yy-\5ldvv] , (5.41) 

X 

with 5l = glr-g^lN, (5.42) 

where in both cases 

75/^ + D75 = aD-i^D . (5.43) 

5.4 The discrete Gross-Neveu model 

Note that for gy = Q and (5p = g^ (15.271) becomes the action of the discrete 
Gross-Neveu model 

5DGN = a'$^ {^(Z} + mo) ^-1/055} • (5.44) 



The model is invariant under finite chiral transformations ip — > ^^75, 'ip 
—75?/' and under a hidden 0(2A^) symmetry. This invariance becomes evident 
in the Majorana representation [76]. Since its beta- function is well known 
in perturbation theory we will use this model as a cross-check of our 
calculation. 



5.5 Schrodinger functional of the CGN model 

The Schrodinger functional of free Wilson and Ginsparg- Wilson fermions in 
two dimensions reads (cf. Section Sj) 

Zo = yDV-D^expl-a^^ V^(a;)D^/'(x)| , (5.45) 

where the Dirac operator D in presence of the boundaries is given by the left 
hand side of (14.351) and (14.48p respectively. The subscript of Zq is to indicate 
that it refers to the free theory. Remember that only the fermion fields at 
times Xq = a, 2a, . . . ,T — a are integrated over in the functional integral and 
that they obey Dirichlet boundary conditions (l4.2H4.3p . It is convenient to 
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assume, as we do in f l5.45p . that the fermion fields are defined at all other 
values of xq as well, but that they are zero there. 

Since the four fermion operators have mass dimension two, they are ir- 
relevant in the discussion of the naturalness of the Schrodinger functional 
boundary conditions (Section [4.1.2l) and the needed boundary counter terms 
(Section [4.51) . Thus with the fermion fields defined as above the Schrodinger 
functional of the chiral Gross-Neveu model is 



DipDip exp 



ij{x)Dij{x) + U2ciOi 



(5.46) 



where the sum is over the four fermion operators multiplied by the corre- 
sponding coupling constant. This notation covers all the possible choices of 
operators presented in the Section 15.31 



5.5.1 Generating functional 

In Section [44l we outlined how to compute expectation values of polynomials 
O in the fermion and anti-fermion bulk and boundary fields. A definition 
of the lattice boundary fields is given by f l4.67p and f l4.68p . The generating 
functional for the interacting theory is then 



Z[r], = J D^'DV' exp < —a^ 



"(5.47) 

Replacing in O the field operators by functional derivatives as defined in 
(I4.7ip . the expectation value can be written in a compact form 

0) = \^OZ[%A . (5.48) 

L ^ J V,r,=0 

For our perturbative expansion it is convenient to rewrite (15.481) . If we also 
replace the fields in the four fermion operators Oj by functional derivatives 
the generating functional becomes 

2 

Z[r],r]]=e^^-'''°' Zo[T],ri]. (5.49) 

where Zq \fj, t]] is the generating functional of the free theory (cf Section I4.4p . 
After integration over the fermion fields it is given by 

Zo[v,v] = Gxp <a^^r]{x) S{x,y)r]{y) \ , (5.50) 

I x,y ) 
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where S{x, y) is the propagator associated with the Dirac operator in (I5.47p 
DS{x,y) = ^5^^y, 0<x,y<T, (5.51) 

CI 

with boundary values 

' 0. (5.52) 



P+Six,y) =P^Six,y) 

xo=0 

Finally we use (I5.49p in f l5.48p to write 



xo=T 



(0) = /loeTE.^.o,^^[-^]j , (5.53) 

L ^ J rj,ri=0 

Expanded in powers of the couplings the generating functional f l5.49p is 
the sum of all vacuum diagrams, that is, diagrams with no external lines. 
Clearly the factor 1/Z in f l5.53p cancels all diagrams in the expansion of the 
expectation value, that contain vacuum diagrams as subdiagrams. Therefore 
we write the expansion of the expectation value as 

(O) = (O)^ + ^J2cj{0 Oj{x))^ + 9tj2cjcj{0 Ojix) Oj{y))^ + . . . , 

X x,y 

(5.54) 

where (^O -^)q, with X a product of four fermion operators or unity, is defined 
as 

{OX\ = l^OXZo[%v]] - 

^ '0 \^Zq J^^_o vacuum subdiagrams 

The terms in the expansion f l5.54p are given by the sum of all Wick con- 
tractions of the fields in the operator O and the increasing number of four 
fermion operator insertions, that do not contain vacuum diagrams as subdi- 
agrams. The basic contractions are given by (I4.73p - fl4.8ip . The contractions 
are most conveniently represented by Feynman diagrams. The rules for the 
corresponding expressions are given in the next section. 



5.5.2 Feynman rules 

Because translation invariance is broken in the time direction the Feynman 
rules in the Schrodinger functional are given in a half Fourier transformed 
space. That is, the half transformed Propagator S{xQ,yQ,pi) is defined 
through 

S{x,y) = ^Y1 ^'''^'''-''^S{xo,yo,Pi)- (5.56) 
pi 
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The rules are listed in Fig. I5.1[ In the case of Wilson fermions the prop- 
agators from the boundary to the interior H{xo,pi), H'{xo,pi) and from 
boundary to boundary K{pi), K'{pi) can be evaluated to some extent (see 
Appendix ETI]). 

Statistical factors Because the four fermion interactions consist of a prod- 
uct of two identical bilinear operators, each insertion of a four fermion inter- 
action can be connected to the rest of a diagram in two ways that lead to the 
same final contribution. In other words there are two different contractions 
giving the same diagram. This factor is always cancelled by the 1/2 that 
comes with each insertion (see f l5.55p ). 

The factor 1/n! from the Taylor expansion needs some more words. At 
the nth order of this expansion it multiplies terms with n four fermion inter- 
actions. Their general strucure is 

(h'o!) = E „,!,JL or... or. (5.57) 

\/=l / mi+m2-\ \-mp=n 

where F is the number of different four fermion operators in the action 
{F = 3 for Wilson, F = 2 for Ginsparg- Wilson fermions). The number 
of equal contractions for a term of this sum due to interchange of vertices 
is mi! . . . mi?!. Thus, taking also into account the overall factor 1/n! 
from the Taylor expansion of the exponential function (see (I5.55P ) we have 
(mi!m2! . . .mF\)/n\ which is exactly cancelled by the factor in Eq. (|5.57p . 
Therefore none of these factors appears in the Feynman rules. 

Sign The sign of a diagram is determined by the number of traces. Each 
trace comes with a minus sign. This gives the overall sign 

/ #traces 
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TflPi ^ a;o H{xo,pi) = S{xo,a,pi)P+ 



^rPi ^ ^° ^Rl5H^xo,pi)j5=riiP_S{a,xo,pi) 

c- " -c 

Ffl Pi p[ Tn' Tr, 6p^p,^K(pi) = Tr, 5p,p, P+S(T - a, a,pi)P+ 



^RPi p[rn' rn^55p,p'^K^ipi)^5 = rnSp^p>^P-Sia,T-a,pi)P- 

TrPi j rRL(j,i)=TnP_S{a,a,pi)P+ 

fci ~ 
xo — — yo J2k^ •5(a;o, yo,ki) 



^]^:><^ E,„.7 f^(^l-^2-fc3+fc4) 

k2 ^ ^ fc3 



• - Fj, for example 75 A" 



i-^>— ^ Fx, current /density insertion, e.g. 7075 



same vertex but separate Dirac and flavour traces 



Figure 5.1: Feynman rules for Schrodinger functional of the chiral Gross-Neveu model. 
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Chapter 6 

Chiral symmetry restoration 



In this Section we use the chiral Ward identity to restore chiral symmetry 
in the Wilson discretisation at finite lattice spacing (up to 0(a)). We show 
here that this can be achieved by a perturbative computation of the critical 
mass rric and a symmetric 5p^s- The value of the third coupling Qy has not 
to be tuned. 

On the lattice, using the Wilson discretisation, chiral symmetry is broken 
and nothing in general prevents mo to take a finite value or (5p 7^ 0. However, 
as we pointed out in Section 13.41 the renormalised axial current (^/x)^ of 
properly defined fields is expect to obey 

{{O)^d,{A,)^{x)) = 0{a). (6.1) 

This condition can be used to fix the bare parameters in perturbation theory 



5ll, lDJ (as well as in numerical simulations (3j|). In the following section we 



discuss the strategy and the result in detail. 



6.1 Correlation functions 

We define correlation functions in the Schodinger functional (SF) set up in 
order to utilise (16. ip . The correlation functions 

fx{x,) = -^J2{ ^(^) H^) CiVi) 75 C(^i) > , (6.2) 
yi,zi 

= 7o75 , Tp = 75 , (6.3) 

to be considered are correlators of a zero momentum pseudo scalar boundary 
state built from the boundary fields (, ( (eq. (I4.67P ) and insertions of the 
time component of the axial current (/^) and the pseudo scalar density (/p) 
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Figure 6.1: Tree level diagrams for fx- 

respectively. The vacuum expectation value (■) has been defined in Section 

As discussed in Section 15.3.11 there is some ambiguity in the operators of 
the four fermion interaction. Here we use the following form of the action 

'S'cGN,w = «^ X] (^w + mo) 

X 

- y\0'ss - O'pp) - \8'lO'pp - ylO'yy] , (6.4) 

for the actual computation. Using (j6.4l) in the perturbative expansion the 
computation gets simplified, because a number of contractions vanish due to 
vanishing fiavour traces. Results obtained for the primed couplings are then 
translated back into the more common unprimed couplings of (15.270 using 
(E29D-(E3l|). 

For small couplings the expectation value in (16.21) can be expanded as 
indicated in (15.541) . Then the right hand side of (16.21) is a sum of an increas- 
ing number of insertions of the interactions but the expectation value taken 
with the free action 5*0 only and all contractions with vacuum subdiagrams 
subtracted (cf. (I5.55P ) 

fxixo) = fPixo) + 5^ c, /i^i(xo) + ^ cjcj /|i,(xo) + 0{c') , (6.5) 

/ i,j 

where 

ci = {g'\5'l,g'l} and Oj = {O'ss ~ O'pp ,0'yy} . (6.6) 
The tree level amplitude 

/f M = E (^(^)rx^(:^)C(l/i)75C(^i) >o, (6.7) 

yi,zi 

is the sum of the two diagrams sketched in Fig. 16. li In these diagrams the 
dotted lines represent the time slices xq = and xo = T. Plain lines are used 
for the fermion propagator and the small open circle in the middle symbolises 
the insertion of a current or density operator. Using the explicit form of the 
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Figure 6.2: First order diagrams for fx- 

porpagator (14.421) the tree level amplitudes can be calculated analytically to 
some extent. The somewhat lengthy expressions are listed in Appendix IB. 31 
The first order amplitudes 

/S(^o) = -^ E {i^i^)^xHx)ayi)i,az,)o'j{u))^, (6.8) 

are sums of five diagrams. Due to the 75-hermiticity of the Wilson-Dirac 
operator, and thus the fermion propagator, the diagrams 1 and 2 in Fig. 
16.21 are equal. In these diagrams the small filled circles (dots) represent the 
insertion of a four fermion interaction. Since the free propagator is diagonal 
in fiavour space the fiavour traces factorise from the rest of the computation. 
Diff'erences in the resulting factors originate from different orders of the A- 
matrices in the four fermion interaction and the number of separate traces. 
Only those combinations of the A-matrices that can be reduced to the identity 
give a non-zero contribution. In the case of the first order diagrams there is 
only one possible order 

A"A" = ^ ^ ' 1. (6.9) 

The number of separate traces is one (diagrams 1,2,5) and two (diagrams 3,4). 
Together with the factor in the definition of the correlation function the 
first order amplitudes can be organised in powers of with p = 0,1 and an 
overall factor (16.91) . 

The second order amplitudes 

/Sj(%) = -^ E (V^Wrx^(x)c(z/i)75C(^i)o;(«)o»>„, 

yi,zi,u,v 

(6.10) 

are sums of the 22 diagrams depicted in Fig. 16.31 The diagrams 1, 5-7, 11 
and 12 have to be counted twice since when refiected at an horizontal line 
they give diagrams with the same numerical value but a different contraction 
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6.2 Strategy and result 
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(just like the first order diagrams 1 and 2 are equal). There are now four 
A-matrices to combine producing factors 

A^A^'A'^A'' 
A"A^ ® A'^A'' 

The number of separate traces reaches from one to three. Together with 
the factor in the definition of the correlation function the second order 
amplitudes can be organised in powers of with p = 0, 1, 2, 3 and an overall 
factor 4(iV2 - l)/N'^. 

6.2 Strategy and result 

We work with the Schrodinger functional of the Gross-Neveu model at fixed 
ratio T/L = 2. The phase 6 = 6i characterising the spatial boundary condi- 
tions (cf. Section 12. 2p is a free parameter of this regularisation. The Ward 
identities are independent of 6. Hence it provides a probe for the critical 
mass in the sense that ( 16.11) must hold for all 6. 

We define the renormalised correlation function (/A(a;o))R and its time 
derivative hA{0,xo/L) as dimensionless quantities 

ifAixo))^ = ZaZI fAixo) and hAi9,Xo/L) = Ldo ifAixo))^^ , (6.14) 

where we introduced normalisation factors for the axial current and the 
boundary fields 

Za = 1 + 0/4'/ + O(c'), = 1 + d Z^'J + 0{c^) . (6.15) 

In general, because (9^y4^(a;) can mix with ^ P(x) we expect for Ha an ex- 
pansion in powers of a/L starting with a linear divergence 

hA{e, xo/L) = A_i{e, xo/L) L/a + Ao{e, xq/L) + 0{a/L) . (6.16) 

Eq. (16.11) enforces the coefficients of the divergence and the finite part to 
vanish 111 

A^i{9,xo/L) = Ao{9,xo/L) =0, for all 9,xo/L, (6.17) 



4(iV2 



1. 



Ar2 

A{N^ - 1; 

iV2 



1 ® 1. 



(6.11) 
(6.12) 
(6.13) 



^Be careful to not confuse the coefficient Ao{6,xo/L) with the time component of the 
axial vector current Aq{x). 
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at amo = amc and 6' p = 6' . (6.18) 

Thus we have a two-dimensional parameter space spanned by 6 and Xq/L for 
which these coefficients must vanish. 

For free fermions (16.11) is satisfied for amo = cf'^c = 0. In the interacting 
theory rric can be expanded 

arric = am^^l cj + am^^jj CjCj + O(c^) . (6.19) 

As we will see also is constrained by (16. ip and is given in terms of the 
other two couplings. 

The dimensionless Ha just as well posses an expansion in powers of the 
couplings 

hA{9, xo/L) = hf + 5^ c, h^}^j + J2 c/cj h^j + 0{c') . (6.20) 

/ I, J 

Now the coefficients of this expansion, sums of lattice diagrams, may be ex- 
panded in powers of a/L. As already indicated /i^'* = ho is at least 0{a^/L^) 
and therefore aml^^ = 0. In the next two sections we compute the divergent 
and finite part of h^^\ and h^^u- Although there will be some analytic ar- 
guments to simplify the expressions, the final evaluation of the contributing 
diagrams is performed numerically as described in Section 12.31 



6.2.1 First order 

Using expansions (16.51) and (I6.19P in (16.20p yields for the first order term 
E ^A,i = E {^1.^ + ^2 + + 2Z^]j) ho) , (6.21) 

with 

ho = LdofA\xo) , hi J = LdofAjM , /12 = q^^o , (6.22) 

all defined at anio = 0. The tree level amplitude ho vanishes identically for 
9 = and is 0(a^/L^) for 9^0. Its derivative /i2 with respect to amo though 
diverges linearly with L/a. This implies at first order 

E c/ {hj + am^^l /12} = 0(a/L) . (6.23) 
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Linear divergence The first order amplitudes hij = Yld=i ^1/ 
sums of the diagrams depicted in Fig. 16. 2[ Only diagrams one to four 
contribute to the linear divergence due to the contact term induced by the 
bubble contraction (cf. Appendix IB. 1.30 . The contribution of diagram five 
hi J is 0(1) and thus is important for enforcing Aq{9,Xo/L) = 0. 

For the linearly divergent diagrams one to four we find (with the help of 



(B.24) 



^^i'^"~2iV ^"^^ ^ 

d=l yi,zi,u 



{ Ao{x) C(2/i) 75 C(^i) H^) {B{uo) + Fi{u) / 9, + ..) ^{u) . (6.24) 

The factor 2(A^^ — 1)/A^ originates from the sum over the A-matrices (cf. Eq. 
(16.91) and comment before). For fixing the critical mass we need only the 
leading term in (I6.24p . Using (|B.21l) and (|B.25l) we find 



d=i 



Thus the contribution of these diagrams is proportional to /12 which also 
multiplies am!-^\ in (I6.2ip . Abbreviating arrrp = am^^] cj the linearly 
divergent part at this order is given by 

A_,ie,xo/L) = (^am« - {s'l + 2g'l^ 2(Ar^- 1) c^,iN,e,xo/L) . 

(6.26) 

Where C^i{N,9,Xq/L) is the coefficients of the linear divergences in /i2 

h2 = C^iiN,e,xo/L)^ + 0{l), (6.27) 

and C^i{N,e,Xo/L) =C^i^o{0,Xo/L) +NC^i^i{e,xo/L) . (6.28) 

The linear divergences in /12 arise from the mass derivative of the two tree 
level diagrams in Fig. 16. 1[ Diagram two receives an additional factor of A^ 
because its evaluation involves two separate fiavour traces (it has two closed 
fermion loops). These amplitudes can be determined numerically. The right 
hand side of (16.261) vanishes for all 9 ,Xo/L only if the critical mass is 



(1) - 1) 

am). ' = — 

A^ 



B, [6'1 + 2g'l^ . (6.29) 
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Finite part Summing up all 0(1) contributions of the five diagrams and 
organising it in powers of N we find that the terms proportional to g'^ and 
g'y cancel. Explicitly we have 

Aoie, xo/L) = ^/^ ^(^^- ^) c,{N, e, xo/L) , (6.30) 

with 

CoiN,e, Xo/L) = Co,oie, Xo/L) + N Co,i{9,Xo/L) . (6.31) 

The right hand side of (I6.30p has to vanish for all 9 , Xo/L. Therefore chiral 
symmetry is restored for 

S% = 0, (6.32) 
at first order of perturbation theory. 

Summary Translating Eqs. (16.290 and (I6.32P back to the unprimed cou- 
plings using (I5.29I) - (I5.31I) yields 

amW = -0.3849001 x {2Ng^ - 4 - 2^y) , (6-33) 
4,. = 0{g') . (6.34) 

At this order we find no constraint on the vector coupling gy. Or phrased 
in another way, the vector coupling has not to be tuned in order to restore 
chiral symmetry. 

6.2.2 Second order 

The relevant terms for the second order term in (|6.20p are 
J2 cicj hf jj = 

c/Cj ^amfjj /i2 + hi J J + am^^Jhsj + aml^Jam^^l + 0(a/L)| , 

i,j 

(6.35) 

with 

KiJ = Ldof'ZjM, h,, = ^h,n h, = --^ho. (6.36) 

It is understood that the first order critical mass is set to the value found 
above and the constraint (5p5 = 0{g^) is employed. Also terms like {Z^\j + 

2Z^^]j) ho and (Z^^j + 2Z^^]) /i^ j that are at least 0{a/L) and suppressed in 
fl635D . 
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Divergences Recall the reduction of bubble diagrams like the diagrams 
1-10 in Fig. 16.31 to insertions of the scalar density and the relation between 
the mass derivative of a expectation value and the insertion of the scalar 
density outlined in Appendix lB.1.3[ After a little thought it is clear that the 
derivatives /is j and cancel all the power divergences in the second order 
amplitudes hijj that are due to bubbles, i.e. the contribution of diagrams 
1-10. That this is indeed the case was checked numerically. The only other 
linearly divergent diagrams are number 11, 12 and 21, 22 in the same Figure. 
In order to keep the equations clear we define 

am^^-* = amfjjCiCj, (6.37) 
IJ 

and the linear divergent part of the summed up diagrams 11, 12, 21 and 22 
E = ^' ^o/^) ^ + 0(1) ' (6-38) 

d=ll,12,21,22 

2 

with CLiiN, e, xo/L) = J] ^" CL{ni0, Xo/L) . (6.39) 

n=0 

The coefficient of the linear divergence at second order is then 

A.,{e, Xo/L) = amf^ C_i + ^^^^^^ Yl ^^^-^ ^-i ' (^-^O) 

i,j 

where C_i = C.i{N,9,xo/L) is defined in (16:271 ). 
Setting the right hand side of f l6.40p zero leads to 



(2) 4(Ar=^ - 1) >^ CLi{N,e,xo/L) 



^ U_,[l\,^,Xo L 



When one evaluates the ratios in the last expression numerically, one finds 
that the 9 and xo/-£/-dependence of the numerator is cancelled by the denom- 
inator. In other words CLUNjOjXq/L) factorises 

CLi{N,e,xo/L) = {D% + D%N) ■ C-,{N,e,xo/L) . (6.42) 

where D^j{ ^ are computable constants. Organising the terms in powers of 
we find at this order 

4(A^2 _ 1) 



- '-^{{D, + ND,){g" + g't.)+2{D, + ND,)g"g'l 

(6.43) 
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with 

1)1 = 0.01195(1), ^2 = 0.22870(1) . (6.44) 

As mentioned the first order result 5'p,s = 0{g'*) must be used here in order 
to cancel all divergences. 

Finite part The first order result allows for a second order contribution. 
Thus we have to include the first order finite part (16.301) at second order 

Aoie, xo/L) = S'p,, ^ilf^ Co + ^^^^7 Yl ^^^-^ ^o" • (6.45) 
The coefficients 

3 

Co"(iV, e, xo/L) = J] AT- C^iie, Xo/L) , (6.46) 

n=0 

receive contribution from all the 21 diagrams of Fig. 16.31 which all have a 
factor of 4(A^^ — 1)/N'^ in common. Cq = Cq{N,6,xq/ L) was introduced in 
f l6.30p and is the finite part of the first order diagrams up to a common factor 
of 2(iV2 -i)/N. Setting flOHll zero yields 

Co{N,e,xo/L) 

Again we find that the 9 and Xo/Iv-dependence of the numerator is cancelled 
by the one of the denominator, which means that chiral symmetry is really 
restored. Note that in the sum on the left hand side many terms vanish 
because of 5'p,s = 0(5'^). Organised in powers of A^ the result is 



6' 



2D. 



P,s 



N 
with 



AT/ + (1 + AT - AT^) + 2g''g% + Oig"') , (6.48) 



L>3 = 0.6192(1) . (6.49) 



Summary Translating Eqs. (|6.43l) and (16.481) back to the more common 
couplings g"^, 6p, gy using ( l5.29p -( !53Tll yields 

am(2) = {D,~N D,) (g' + g^r) + 2{D, - A^ D,) g' gl + 0(/) , (6.50) 
<5p,, = D3 [Ng'-2g^gl-g'y]+0{g^). (6.51) 
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6.3 Conclusion 

We demanded the chiral Ward identity (16. ip to hold up to 0(a) on the lattice 
for the renormalised operators. The computation is carried out in second 
order perturbation theory in the Schrodinger functional of the Gross-Neveu 
model with Wilson fermions. The result is that the bare mass has to diverge 
in order to cancel a linear divergence and that the chiral symmetry breaking 
coupling 6p, although zero at first order, has to take a finite value at second 
order. 

Explicitly the result is 

arric = am^ + am^^^ + 0{g^) . (6.52) 

with 

am(^) = -Bi X {2Ng^ - - 2^y) , (6.53) 

am(2) = {D^ -ND,) (/ + g^.) + 2(^2 - N D,) g^ gl + 0(/) , (6.54) 

and 

5l^, = D, {Ng'-2g^gl-gtr)+0{g'). (6.55) 

We find no constraint on the vector coupling gy. Or phrased in another way, 
the vector coupling has not to be tuned in order to restore chiral symmetry. 



In order to compare with the large result of [SOfl we set and rescale 
g^=gl/N, 5l = {gl + gl)/N and gl^gl/N, (6.56) 
and take the N oo limit in (16.521) and (I6.55P 



am 



Tv^oo _ 0.7698002 x + 0(/) , (6.57) 



[gl/glYr^ = 1 - 0.6192(1) xgl + Oig') . (6.58) 

For amc this is the complete large result. The authors neglect the vector 
four fermion interaction from the beginning, but their result is not changed 
if it is taken into account. In the large N limit chiral symmetry is restored 
for l/gp = l/gg + 0.619 which is reproduced by our result at second order. 
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Chapter 7 

Renormalised coupling 



If the time extension T of the Schrodinger functional of the massless chiral 
Gross-Neveu model is fixed to a multiple of the spatial extension L, say 
T = 2L, then all dimensionfull quantities will depend on L. Keeping the 
volume finite usually leads to systematic errors. But here the finite volume 
is utilised to probe the theory. That is renormalised quantities are defined 
at the scale fi = 1/L. Therefore by definition there are no finite size effects 
in the Schrodinger functional. 

In view of a Monte-Carlo simulation of the model we aim here at the 
definition of renormalised couplings at zero renormalised mass in terms of 
renormalised correlation functions of the boundary fields f l4.67p . (14.681) . For 
each coupling g] in the action a combination Fj of such correlation functions 
defines the corresponding renormalised coupling gj such that it is equal to 
the bare coupling at leading order of perturbation theory up to corrections 
of order a 

~gj = Fj such that g] = gj + Oiag^) + Oig"") . (7.1) 

In the continuum limit of asymptotically free theories this ensures gj "'=' 
gj + 0{g^)^ In turn, as we have seen in Section [2.1.31 such a definition of 
the renormalised coupling ensures the universality of the first two coefficients 
of the corresponding beta-function in the single coupling case. For multiple 
couplings the situation is more complicated (see Section 12.1.41) . Nevertheless 
condition (I7.ip is a possible choice also in that case. 

In the next Section the used correlation functions are discussed. In Sec- 
tion [7i2] they are formally expanded in powers of the couplings exploiting the 
cancellation of some diagrams. The results for Wilson and Ginsparg- Wilson 

^Throughout the thesis 0((g^)") means all terms 0((.gjj"i . . . (5/^ )"""), X^Zli '^^ = 
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fermions are presented in Sections [7. 31 and FTH respectively. In the last Section 
(17. 5p we use the results to compute the ratio of the A-parameters. 

7.1 Correlation functions 

The two- and four-point correlators of the boundary fields are defined as 

/2 = -]^E(^MC'(^i)>' (7-2) 

UlZl 

and 

f^ = -WN^^T)r^ E (c'(^i)75Ar(^i)c(i/i)75A'^c(^i)>- (7.3) 

uiviyizi 

We sum over the spatial direction to project on to zero momentum. The 
powers of the lattice spacing a are chosen to render the correlation func- 
tions dimensionless. The remaining factors are to normalise the tree level 
amplitudes. 

The matrices Tb and contract the Dirac indices of the boundary fields. 
In two dimensions such a matrix can be expanded in the identity, 7^ and 
75. But only terms with Tb and equal to 71 and/or 75 have a non- 
vanishing contribution because of the projectors P± at the boundary of the 
Schrodinger functional (apparent for example in the basic Wick contractions 
at the end of Section or the Feynman rules in Fig. 15. ip . For the same 
reason the contribution of the possible combinations of the two are all equal 
up to a factor of i and the overall sign (^±75^:^ = ±iP±'yiP^). Therefore 
in the computation of we use for the Dirac structure on the boundaries 
T's = Tb= 75- 

For the flavour structure we choose flavour vector bilinears. In combina- 
tion with four fermion interactions of the flavour scalar type (|5.7I) this choice 
simplifles the perturbative computation. This is because diagrams which do 
not connect the two boundaries involve the flavour trace tij {A"} and hence 
vanish. 

The renormalised two- and four-point correlation functions are 

{f,)^ = Zlh and {U\ = ZlU (7.4) 

The normalisation factor for the boundary flelds was already introduced in 
fl6T5|) in Section El 
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In the ratio of appropriate powers of the renormalised quantities the un- 
known wave function normalisation is cancelled. Then it is enough to evalu- 
ate the unrenormalised correlators at zero renormalised mass 

R{e) = - 1 = - 1 at mR = . (7.5) 



(/2)r (/: 



2 



For later convenience we subtract the free theory value at zero mass. As 
indicated the ratio R = R{9) depends on the phase 9 = 9i parametrising the 
spatial boundary conditions. This dependence will allow us to find combina- 
tions R{9) + b ■ R{9') to define the renormalised couplings as in (17. ip . 



7.2 Perturbative expansion 

In order to expand Rg in powers of the couplings gj we need the expansion 
of and /2. We abbreviate c/ = (^f and indicate that opposite to f l6.6p here 
c/ refers to g^, 6p and gy. 

h = + ci ffj + cjcj ffjj + 0(c3) , (7.6) 

and (/a)^ 

(hf = (/f )' + ^- + c,c, {2Sf^Sfh + + 0(c^) . (7.7) 

Multiplying (17.61) with the inverse of (17.71) yields 

/4 _ 1 , J (1) _ ,(1) ,{0) 1 

" (/f)n ' ' J 

+ -rB^- f'S., - VSjfP + 3/l','/£ - 2^ + 0(c») . (7.8) 

(/<°>) I h J 

The first and second order diagrams of /a and /4 are listed in Figs. 17.11 and 
17. 2[ Half of the diagrams of /4 at first and second order can be shown to 
be identical to products of /a diagrams. To see this, consider a diagram /^"^ 
at order n in the expansion of f^, where the subscript i labels the diagram. 
Then we define: 

Definition 1 If a diagram f^^^ at order n in the expansion of fi can he cut 
horizontally without cutting through a fermion line or a vertex, then it is 
called reducible. 
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f(i) _ 

J2 — 



J4 — 



Figure 7.1: First order diagrams of /2 and /4. The bubble diagrams stand for the sum 
of the connected (one trace) disconnected (two traces) diagram. 



f(2) _ 
■12 — 

f(2) _ 
/4 — 





+ 2 




+ 2 




+ 






Figure 7.2: Second order diagrams of /2 and f^. The bubble diagrams stand for the sum 
of the connected (one trace) disconnected (two traces) diagram. 



In Appendix IB.2.41 we prove the following Lemma: 

Lemma 1 If a diagram /I"'' at order n in the expansion of f^ is reducible, 

then f^^j^ can be written as a product f!^^ ■ f^l of two diagrams appearing 

in the expansion of f2, with n = r + s. The amputated part of and 

f2k is equal to the amputated part of the upper and the lower half of f^^^', 
respectively. 

Diagrammatically this is depicted in Fig. 17.31 As a consequence all prod- 
ucts of /2 diagrams in the expansion (17.81) are cancelled by corresponding f^ 
diagrams. Only /4 diagrams that are not reducible are left 

/4 ^ 1 , ^ f (1) 
2 „(0) •'4,/,a 



C£CJ|.(2) 2f(2) I f(2) , f(2) ,4. 

^ (.(0) ] J'i,IJ,a^ '^J4,IJ,h^ J4:,IJ,c^ J4,IJ,d^ ^J^ 



Cj J (2) 

(n^ S .1/1 T T., 



From this result one easily reads off the expansion of R 

R{e) = R^^\e) + R^^\e) + o{c^) , (7.10) 





Figure 7.3: Diagrammatical representation of Lemma [TJ 



7.3 Wilson fermions 
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with 

R^^\e) = ciRf{e) and R^^\9) = cjcj Rf){e) . (7.11) 
The first and second order terms are 

4'^W = ^/Sa, (7.12) 

_ ^/ f(2) , 9f(2) , .(2) .(2) , ..(2) 1 

^IJ — (0) S J4,IJ,a ^i4,/J,b i4,/J,c i4,/J,d 4,IJ,e ^ (o) ^4,J,a f • 

74 t J2 ) 

(7.13) 

We now proceed with the explicit computation for Wilson fermions. 

7.3 Wilson fermions 

The lattice action of the chiral Gross-Neveu model with Wilson fermions we 
use here is given by (15.271) 

'S'cGN,W = «^ X] ] ^ + "^o) ■ip 

X ^ 

- -,g'iOss - Opp) - \5lOpp - \glOvv^ ■ (7.14) 

The renormalised mass vanishes if the bare mass is set to its critical value 
and the coupling of the pseudo-scalar interaction to its symmetric value 

mo = rric and 5p = . (7-15) 

Since Lemma [U holds for any mo, it is also true for mo = rric- In order to 
incorporate the expansion of rric in the expansion of R, two new terms have 
to be added at second order, corresponding to the derivative with respect to 
amo of the first order term in (17.91) {dm = d/damo) 

R?iO) = ^fi]l, (7.16) 



Rf]{e) 



f(2) _|_9fl^i _|_ f(^) I f(^) j_Ai'^^' 
Ji,IJ,a ' '^Ji,IJ,b ' J4,IJ,c ' JAJJA "i" ^JAJJ. 



c(2) 



(2) 



K2) 



(0) 



(1) 



/ff + am\\'j dmf2 
AO) 

J2 



(1); 



c(0) 



(1) 
4,J,a 



(7.17) 
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where all diagrams are evaluated at mo = 0. We do not set 6p = from 
the beginning, but keep it as a free parameter and only set it to its symmetric 
value at the end of the computation. This will allow us to compare the general 
result with computations in the discrete Gross-Neveu model (cf. Section EH]). 

In the case we have to evaluate only one diagram which involves 

a sum over the time coordinate of the four fermion interaction. This sum 
can be computed analytically to some extent (cf. Appendix IB. 2. 21) . i.e. the 
continuum limit can be extracted. Using (IB. 351 1 and ( IB. 271 1 we find 

i?(^)(0) = {2/ + Si {A{e) -l)-2gl B{e)] + 0(a) . (7.18) 



with 



L 



A{e) = ^ sinh(2^r/L) n° 1 (7.19) 

B{0) = cosh(2^T/L) 1 (7.20) 
C{e) = cosh\eT/L) 1 . (7.21) 

(2) 

In the case of Ryj{9) there are five non-reducible diagrams involving 
one momentum loop. The diagrams are evaluated numerically for a range 
of lattice sizes and several 9 values and extract the finite and logarithmic 
divergent terms (cf. Section [231 and Appendix IB.2.3p . Using Eqs. (|B.39l - 
[Ril we find 



R^'\0) = ^^1^ { - 4^^ (AT - Bi9)) + 24 (AT - 1) iAi9) - 1) 

-Ag^Sl iiN-l)Ai9) + Bi9)~N) + 8g'ySlAi9)] + ... , (7.22) 
where the dots indicate the finite part that has been suppressed here. 



7.3.1 Discrete Gross-Neveu model 

Our general four fermion interaction theory also contains the well studied 
Gross-Neveu model [I^. We call it here the discrete Gross-Neveu model for 
its discrete chiral symmetry ip TsV'? ^ ~^ and in order to discriminate 

it from the chiral Gross-Neveu model. The discrete Gross-Neveu model has 
only the scalar four fermion interaction g'^ /2 (ipip)"^ (cf. Section [5^ . In our 
notation this amounts to setting 



(7.23) 



7.3 Wilson fermions 
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This model possesses a 0(2A^) symmetry that allows no other interaction 
term. This strictly holds only for 6^ = , vr. For all other values the bound- 
ary conditions break this symmetry. However, the local Ward identities as- 
sociated with this symmetry will still hold and therefore the ultra-violet 
divergences are expected to remain unchanged. 

The results of the last section can be used to calculate the one-loop beta- 
function, which then can be compared to previous results. 

With the settings (E23D, eqs. (TTlSll and (17:2211 simplify to 

^S'^^^^ = 2LW) ^^^^^ ^ ^ ^ ^^""^^ ' ^^'^^^ 

< W = ^ (N - 1) (Aie) + 1) + . . . . (7.25) 

A renormalised coupling can be defined by 

algn = Fdgn = rpl^^^-^^^l^ RdgniO) , (7.26) 

with the expansion 

4„ =° 9' - 9' ln(a/L) - Cdgn} + 0(/) , (7.27) 

with Cdgn = Cdgn + c^gjj N. The corresponding beta- function (12.71) is (/i = 1/L) 

This result is in accordance with previous continuum calculations 
78 and recent lattice calculations 76 



7.3.2 Chiral Gross-Neveu model 

Now that we have confidence in our expansion, we set the pseudo-scalar cou- 
pling to its symmetric value 6p = 6p ^ to ensure continuous chiral symmetry 
at finite lattice spacing (up to 0(a), cf. Chapter [6]) 

Sl^, = D, {Ng' -2g^gl- gt.) + 0{g') . (7.29) 

Since the correction to zero is of order g'^, this enters as a finite term in 
R^'^\9). Explicitly, in the chirally symmetric case Eqs. (j7.18l) and (I7.22p 
become 

^SW = ^^ {/-^7yi?W} + 0(a), (7.30) 
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RSM = -Y§^^ (N - B{e)) g^ + ..., (7.31) 

where we suppressed the finite part. Thus renormalised couplings may be 
defined through 

f^Fg= ^^^^^^ _ {B{e) i?egn(0) - C{e) i?egn(^)) , (7.32) 

and 

fv^Fv = ^^^^^^ _ (/?cgn(0) - C{e) /?eg„(^)) . (7.33) 
They obey the condition (17. ip as can be seen from their expansions 

f g^-g^^ ln(a/L) + c, + 0(^7^) , (7.34) 

TT 

and 

fv =° al-g'- ln(a/L) + cy + 0(/) , (7.35) 

TT 

where the finite parts at second order are given by 

4 I 4 I 2 2 D^D(6) /,^4 r)22 4\ 

C9 = Cg,ggg + Cgyv gy + Cgyg gyQ ^ — [N g -2g gv-9v) > 

(7.36) 

and 

4 I 4 I 2 2 D-^D{6) /,^4 r)22 4\ 

cy = cv,gg g + cy,yy S-y + Cvyg 9v9 ^ (^5' -^9 9v ' 9v) > 

(7.37) 

where D{9) is the ratio 

The coefficients cj^jk = c^i jk^ ^} jk ^ show a rather mild dependence on 6. 
They have been determined numerically for several values of 6 in the interval 
[0, 1]. The result is shown in Fig. 17.41 and Fig. 17.51 where open and filled 
symbols refer to and c] jj^ respectively. The last term in (17.361) and 

( I7.37P originates from 6p^s (I7.29p . which has to be included at this order. 
It is straightforward to derive the associated beta- functions (/x = ^/L) 

Pg{f) = ~f- + Oif), (7.39) 

TT 

and 

Pvi~9') = ~~9' - + Oif) . (7.40) 
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Figure 7.4: Dependence of the finite part of on 9 at order 
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Figure 7.5: Dependence of the finite part of Qy on 6 at order g'^. 
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Thus, from f l7.39p we see that the coupling is asymptotically free and from 
f l7.35p it is obvious that the coupling gy receives an additive renormalisation. 



The one-loop beta-functions derived here agree with the ones derived in the 
MS scheme in Ref. 



22 



As we have seen in Section EH] there is some freedom in the choice of four 
fermion interactions due to the identities f l5.1ip -f l5T3ll . In particular, it is 
possible to find a combination of terms were one of the two beta-functions 
vanishes. This is also indicated by the formal continuum argumentation in 



2A l75|. Indeed if we use flBTal) to rewrite fTTTiD as in (02D . the new 



coupling 5y is related to the original ones by (15.331) 

5l = gl-gVN. (7.41) 

Then fl7.32p remains unchanged. But the renormalised coupling associated 
with 5y is given by 

~5l=fv- f/N = 6l + cy- cjN + 0{g') . (7.42) 

There is no logarithmic divergence and therefore the corresponding beta- 
function vanishes at this order 

Ps{K) = OCg') . (7.43) 



7.4 Ginsparg- Wilson fermions 

The free modified Neuberger Dirac operator defined in Section [1.3.1l obeys 

75Dn + ^n75 = a^N75^N + As, (7.44) 

The parameters c and s are set to 1 and 0, respectively, in the following. The 
term is supported in the vicinity of the boundaries and decays exponen- 
tionally with the distance to them (cf. Fig 14. II) . In particular, the rate of the 
decay is constant if the distance is measured in lattice units. Then transfor- 
mation (I5.34P is a symmetry of the action and the associated Ward identities 
are expected to hold in the interior (well separated from the boundaries) of 
lattice with small corrections which vanish in the continuum limit. Therefore 
we can use an action with two four fermion interactions (cf. Section I5.3.2p 

ScG^,Gw = a'Y^ D^i' - b'iOss - Opp) ~ Ig^yOvv} , (7.45) 

X 

to compute the correlation functions (17. 2p and (17. 3p in the Schrodinger func- 
tional with Ginsparg- Wilson fermions. The operators Oj differ from the 
operators Oj by the substition ip ^ ip = {1 — ^D)ip. 



7.4 Ginsparg- Wilson fermions 
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Figure 7.6: Cut-off dependence of the coefficient of the leading order term in the ex- 
pansion of the renormaUsed coupUng of the discrete Gross-Neveu model. It is unity in the 
continuum Hmit, as it should be. For comparison we plot Ginsparg- Wilson and Wilson 
fermions. 



The free propagator entering in the perturbative expansion is computed 
numerically as indicated in Section 14.3.21 for lattice sizes L/a = 4,5,... ,48 
and e = 0,0.1,0.5,1. 

7.4.1 Discrete Gross-Neveu model 

As in the case of Wilson fermions we check if the beta-function of the discrete 
Gross-Neveu model is correctly reproduced. The renormalised coupling was 
defined in (j7.26l) . At finite lattice spacing one expects 

4n = 9^ fcdgn(a/L) - / ln(a/L) - Cdgn(a/L)| + 0(/) , (7.46) 

where CAgn{.o,/L) = c'^^^^a/L) + c^^^i'^/L) N. The coefficient of the leading 
order term kdgn{a-/L) is expected to be unity in the continuum limit. In Fig. 
I7.6l its lattice spacing dependence is plotted for 6* = , 0.1 , 0.5 , 1. In any case 
it has the right continuum limit (the systematic error of the extrapolation is 
O(10~^)). In Fig 17.71 the finite part at second order is plotted for three 9- 
values. These numbers are obtained by subtracting the expected logarithmic 
divergence from the second order diagrams (we also determined the coefficient 
of the logarithmic divergence directly with results compatible with — (A^ — 
l)/7r and a systematic error of O(10~^)). Note that the cut-off effects of 
the Ginsparg- Wilson fermions are roughly twice as large as for the Wilson 
fermions. 
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Figure 7.7: Subtracting the expected logarithmic divergence (as indicated in l|7.46p ) 
from the second order diagrams one obtains the finite part of g\^^ at order g"*. We plot 
here the result for three different values of 6. Open and filled symbols refer to c^^^{a/L) 
and (a/L) respectively. 



Thus the coupling of the lattice theory with the modified Neuberger op- 
erator (j4.48l) can be renormalised in the same way as with Wilson fermions. 
We now proceed with the chiral Gross-Neveu model. 



7.4.2 Chiral Gross-Neveu model 



The renormalised couplings for the two interaction terms are defined in the 
same way as in the Wilson case, i.e. by implementing the ^-dependence of 
the leading order term in the expansion of R{0) 



f^F,= ^^^^^^ _ {B{d) R„{0) - C{e) R„{e)) , (7.47) 
and 

9v = Fv = y(-2?(6') - 1) ^^"^""^^""^^^ ~ ^^^^ -Rcgn,gw(6')) • (7.48) 

The difference is that no additive mass renormalisation is needed (amo = 0) 
and that the pseudo-scalar coupling vanishes exactly ((5p = 0). 

Since in the case of the modified Neuberger operator we have no analytic 
handle on the leading order, we first check whether the definitions above 
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Figure 7.8: Cut-off dependence of the coefficients of the leading order terms in the 
expansion of the renormalised coupling of the chiral Gross-Neveu model for 9 = 0.5. 
For comparison we plot Ginsparg- Wilson and Wilson fermions. 

satisfy (17. ip . To tliis end we expand 

f = / kg^g{a/L) + gl kgy{a/L) + 0{g^) , (7.49) 
and 

= kv,g{a/L) + gl kv,v{a/L) + 0(/) . (7.50) 

The coefficients were computed for lattice sizes L/a = 4,5,... ,48. Fig. 
17.81 shows the coefficients of g"^ and Fig. 17.91 the coefficients of gy. Their 
extrapolations to the continuum limit have the expected values within the 
systematic errors. For example, in the case 9 = 0.5 the result is 

kg^g = 0.999995(17) kgy = 0.000012(27) , (7.51) 
and 

kyg = -0.000005(18) kv,v = 1.000011(28) . (7.52) 

The results are similar for = 0.1 , 1. Again the cut-off effects of Ginsparg- 
Wilson fermions exceed the ones of Wilson fermions. 

At next to leading order we expect to find the logarithmic divergences 
with coefficients —N/n and — l/vr, respectively, for g^ and gy. As in the 
case of the discrete Gross-Neveu model we first determined the coefficient 
of the divergence explicitly. After we were convinced that it has the right 
value, we subtracted it from the sum of second order diagrams. The resulting 
amplitudes should have a finite or vanishing continuum limit. Expanding the 
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Figure 7.9: Cut-off dependence of the coefficients of tlie leading order terms in tlie 
expansion of tfie renormalised coupling gy of the chiral Gross-Neveu model for = 0.5. 
For comparison we plot Ginsparg- Wilson and Wilson fermions. 



renormalised couplings to second order 

f = / kgML) + gl kgy{a/L) - ^ \n{a/L) + Cg{a/L) + 0(/) , 

(7.53) 

and 

fv = / kv,g{a/L) + gl kv,v{a/L) - / - ln(a/L) + cy(a/L) + 0(/) , 

TT 

(7.54) 

these amplitudes are Cg{a/L) and cv{a/L). They can be sorted by powers of 
the bare couplings 

= Cg^gg 9^ + ^gyv Qv + Cgyg gW > (7-55) 
and 

cv = cv,gg g^ + cvyv Qv + ^vyg Qvg^ , (7.56) 

and powers of N 

ci,jK = cIjk + c] jK N . (7.57) 

In the continuum limit we do not expect to find the same values for these 
coefficients as in the Wilson case, because the finite part of the renormalised 
coupling is scheme dependent. However, as indicated above they should have 
a finite or vanishing continuum limit. Any uncancelled divergence would 
mean that additional terms are needed in the action to absorb them. But 
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Figure 7.10: Cut-ofF dependence of the coefficients of the second order terms in the 
expansion of the renormalised coupling g'^ of the chiral Gross-Neveu model for 9 = 0.1. 
Open and filled symbols refer to (Pj jj^{a/L) and c\ jj^iajV) respectively. 
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Figure 7.11: Cut-off dependence of the coefficients of the second order terms in the 
expansion of the renormalised coupling gy of the chiral Gross-Neveu model for 6 = 0.1. 
Open and filled symbols refer to Cj jj^{a/L) and c\ jj^{a/L) respectively. 
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Figure 7.12: Dependence of the continuum extrapolation of the finite part of 5^ on 9 at 
order g'^. The symbols are Ginsparg- Wilson and curves are Wilson fermions. Open and 
filled symbols refer to Cj and c\ jj^ respectively. 



as can be seen from Figs. 17.101 and 17.111 there is no such divergence, all 
coefficients have a well defined continuum limit. The dependence of the 
values of the extrapolation on 6 is depicted in Figs. I7.12l and l7.13i Although 
the amplitudes are different, the functional dependence is, as it should be, 
identical to the Wilson case. 

Since we find the same divergences the beta-functions associated with the 
renormalised couplings are the same as in the Wilson case. This, of course, 
comes as no suprise. If the theory with the modified Neuberger operator 
in the action is renormalisable and describes the same continuum theory, all 
coefficients of the beta-function must be equal since it is universal as a whole. 



7.5 Summary and ratio of Lambda parameters 

In the last two Sections we calculated the renormalised couplings and beta- 
functions of the discrete and chiral Gross-Neveu model in one-loop pertur- 
bation theory for Wilson and Ginsparg- Wilson fermions in the Schrodinger 
functional. Both models have a asymptotically free coupling. In the chiral 
Gross-Neveu model there is also a coupling that does not renormalise, i.e. its 
beta-function vanishes. The first coefficient of the beta-function is universal 
and we correctly reproduce the known values 

N -I N 
DGN: &o = CGN: 60 = • (7.58) 

7r vr 
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Figure 7.13: Dependence of the continuum extrapolation of the finite part oi (jy on 9 at 
order g"^. The symbols are Ginsparg- Wilson and curves are Wilson fermions. Open and 
filled symbols refer to Cj jj^ and c} respectively. 

A nontrivial test of the equality of the theory with Wilson and with Ginsparg- 
Wilson fermions is now the computation of the ratio of the lattice A-param- 
eters to show its independence of the angle 9. 



- 




□ 
o 

V 


cvyv 

Cv,gV - 

Cvyv ' 

Cv,gV 




V 




V 








-G 



0.2 0.4 0.6 0.8 1 

e 



Discrete Gross-Neveu model: The continuum limit of the renormalised 
coupling of the discrete Gross-Neveu model calculated with Wilson fermions 

f =° 9l + 4 {bo Ha/L) + c-} + Oigl) , (7.59) 
and with Ginsparg- Wilson fermions 

g' =° + 9^ {bo Hc^/L) + + 0{gl) , (7.60) 

must be equal. This allows to relate the bare couplings, e.g. 

9l^ = 9l + aigt + 0{gl), a, = - c^^ . (7.61) 

The finite parts c^^° and Cg^° are functions of 9. Since it is a probe to 
the theory like external momenta, the relation between the bare parameters 
(17.6ip can not depend on it. Therefore the difference ai has to be independent 
of 9. This is equal to saying that the ratio of the lattice A-parameters 

Alat, gw/ Alat, ^ = { ^ } • C^-^^) 
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has to be independent of 6. The lattice A-parameter is analogous to the 
A-parameter in f l2.24p . but with the renormalised coupling replaced with the 
bare coupling and the beta-function replaced with the lattice beta-function. 
For the ratio of the lattice A-parameters in the discrete Gross-Neveu model 
at ^ = 0, 0.5 and = 2, 4 we find 



A 



LAT, gw/ AlAT, w 



N=2, 



=0 [0.5] 



A 



LAT, gw w 



0.987(7) [0.986(12)] 



0.712(12) [0.713(21)] 



N=4, 9=0 [0.5] 



(7.63) 
(7.64) 



Chiral Gross-Neveu model: The chiral Gross-Neveu model has two cou- 
pling constants and the result for the renormalised couplings can be given in 
the following form for Wilson fermions 



~2 a-^0 2 

9 - " 



N 



gi-gl-\n{a/L)+c^ + 0{g: 

71 ^ 



?2 e-2 



and Ginsparg-Wislon fermions 



?2 a^O (-2 



K,,. + cV-cr/N + 0{g: 



gwJ 



(7.65) 
(7.66) 

(7.67) 
(7.68) 



These are couplings of (15.241) . The finite parts are composed of three terms 



~ "^9,99 9i + '^g,VV 9v,i + '^l,Vg 9v,i9i ; 
Cy = Cy^gg 9i + Cyyv 9v,i + '^V,Vg 9v,i9i 



(7.69) 
(7.70) 



for z = w, gw (no sum). The renormalised quantities must coincide. This 
yields relations between the bare couplings 



^Jw = 91 + ^9 + 0{gl) , a 
4,w + «y + 0«J, 



ay 



ag/N. 



(7.71) 
(7.72) 



Again the coefficients Qg and ay have to be independent of 9. That this is the 
case can be seen in Fig. 17.121 and 17.131 There and ay are the differences 
between the curves and the corresponding symbols. It is finite, but constant. 

The beta-function of the vector-vector coupling 6v vanishes at one-loop 
(I7.43p . This is in agreement with two-loop calculations [22| and formal con- 



tinuum arguments |74l . |75|. Furthermore the vector- vector coupling was 



7.5 Summary and ratio of Lambda parameters 



89 



not constrained by the restoration of chiral symmetry in the case of Wil- 
son fermions (see Chapter [6l). Therefore no value of the bare coupling Sy 
is distinguished. For each fixed value the theory is an one coupling theory. 
The ratio of the lattice A-parameters can then be calculated as in the dis- 
crete Gross-Neveu model (Eqs. (I7.6ip and (17.621) ). For convenience we take 
^y,w — 0- From (17.721) it follows that 5y^^^ = 0{g^) and the coefficient ag 
simplifies to 



'g,gg 



- C 



.gw 
'9,99 ■ 



(7.73) 



We list the ratio for 9 = 0.1, 0.5 



ALAT,gw/ A-LAT,' 



,w 



0.9893(7) [0.9892(10)] 



(7.74) 



Ar=2,6»=0.1 [0.5] 



/ A-LAT,' 



0.776(1) [0.776(2)] . 



(7.75) 



Ar=4,6»=0.1 [0.5] 



The ratio of the A-parameters has been calculated in QCD on a lattice with 
periodic boundary conditions in j79|. The values there are similar to the ones 
obtained here. 
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Chapter 8 
Conclusions 



The main result of this work is the formulation of the Schrodinger functional 
(SF) for fermionic models of the Gross-Neveu type with a finite number N 
of fermion flavours. In 1-loop lattice perturbation theory we showed that the 
theory is renormalisable with Wilson and Ginsparg- Wilson fermions. This is 
the flrst check of the recently proposed Dirac operator (32| beyond the free 
theory. 

In two dimensions four fermion interactions have dimensionless coupling 
constants. We discussed the symmetry properties of the four fermion interac- 
tion terms and the relations among them. Due to its Abelian chiral symmetry 
the continuum chiral Gross-Neveu (CGN) model has two independent inter- 
action terms. A possible choice is {{ipip)^ — (^'75^')^) + 5*^/2 ii^lfii^)^- 
Because Wilson fermions explicitly break chiral symmetry, the most gen- 
eral lattice action for the chiral Gross-Neveu model has an independent four 
fermion interaction that breaks chiral invariance, for example 6p/2 (t/j'y^i/j), 
and a mass term. The model with three couplings becomes the discrete 
Gross-Neveu model (DGN), when two of the couplings are set appropriately. 

Like in QCD chiral Ward identities can be used to deflne the critical 
mass. The Ward identity of the local axial current, strictly valid only in the 
continuum, is demanded to hold on the lattice up to corrections of 0(a). The 
critical mass cancels the linear divergence that appears in the Ward identity 
because of operator mixing. However, for the restoration of chiral symmetry 
it is also necessary to tune the symmetry breaking coupling to its symmetry 
restoring value. We computed the critical mass and the symmetric coupling 
up to second order in lattice perturbation theory. We find no constraint on 
the vector- vector coupling gy. Or phrased in another way, the vector- vector 
coupling has not to be tuned in order to restore chiral symmetry. This is the 
first determination of this parameters at finite- A^. The result is consistent 
with calculations in the large- A^ limit. 
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Renormalised couplings for and gy are defined at vanishing renor- 
malised mass. The natural choice are boundary-to-boundary correlation 
functions with four external fermions. For Wilson fermions the non-trivial 
dependence of the first order diagrams on the phase 6, which parametrises 
the spatial boundary conditions, could be calculated analytically. This was 
crucial in order to define the renormalised couplings in such a way, that 
they are equal to the corresponding bare coupling at leading order pertur- 
bation theory. This definition has then been applied in the calculation with 
Ginsparg- Wilson fermions, where one has no analytic handle on the diagrams. 

The coefficient of the logarithmic divergence and the finite part at 1-loop 
is computed numerically to a high precision using the method described in 
(45I . We correctly reproduced the first coefficient of the beta- function for the 
CGN and DGN. Both models have a asymptotically free coupling. Further- 
more the couplings can be redefined such, that the vector-vector coupling 
does not renormalise, i.e. its beta-function vanishes. This is consistent with 
2-loop calculations and formal continuum arguments. The finite part shows 
a mild dependence on 6. The cut-off effects are clearly 0(a). 

The definition of the renormalised couplings is suitable for a computation 
in Monte Carlo simulations. The boundary-to-boundary correlation func- 
tions are easily implemented and the 6 dependence is strong enough to dis- 
criminate the couplings. Since the SF is a finite size regularisation scheme 
the spatial extension provides a natural scale in the system and it is possi- 
ble to define step scaling functions jl| for the renormalised couplings. The 
renormalisation group invariant step scaling functions can then be used as 
benchmark observables for universality studies of different lattice actions. 

We used the results of the computation with Wilson fermions to study 
a recently proposed modified Neuberger-Dirac operator |32| in 1-loop lattice 
perturbation theory. The operator is compatible with the SF boundary con- 
ditions in the sense of Section [4.1.2[ The operator of the free theory is shown 
to be local (but not ultralocal) and to obey the Ginsparg- Wilson relation 
up to terms localised at the boundaries with exponentially decreasing tails. 
Thus the lattice chiral symmetry associated with Ginsparg- Wilson fermions 
is a symmetry in the interior of the lattice and correlations of local fields 
at physical distances from the boundaries obey the same Ward identities as 
they do on periodic lattices. After the substitution ip ^ {1 — ^D)ip the new 
four fermion operators transform under the lattice chiral symmetry in the 
same way as the old ones under the continuum symmetry. Therefore 6p = 
is protected by the symmetry and no tuning is necessary. 

The modified Neuberger-Dirac operator (14.481) is rather complicated. The 
eigenfunctions of the operator under the square root were derived. We did 
not succeed in deriving an analytic expression for the Dirac operator. But the 
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Dirac operator and its inverse can be computed numerically using standard 
techniques [72 . 

The computation of the previously defined renormalised couplings was 
repeated with Ginsparg- Wilson fermions. For the coefficients of the loga- 
rithmic divergence we found perfect agreement with the Wilson result, thus 
proving the SF with Ginsparg- Wilson fermions to be renormalisable at the 
one loop level. As expected, the regularisation dependent finite part differs. 
The cut-off effects are 0(a) and larger than in the Wilson case. To achieve 



0(a)-improvement one would have to redefine the correlation functions [80 
and to introduce four fermion interaction terms at the boundaries. 

Finally we used the 1-loop calculation of the renormalised coupling with 
Wilson and Ginsparg- Wilson fermions to compute the ratio of corresponding 
A-parameters, which yields reasonable results. 

With this work the Schrodinger functional for Gross-Neveu models is 
well established. It can be used as a benchmark system for fermion actions. 
It should be possible to simulate with Wilson as well as Ginsparg- Wilson 
fermions. 
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Appendix A 
Notation 



A.l Definitions 
A. 1.1 Dirac matrices 

The 7-matrices are defined through the Clifford algebra in Euclidean space 

Since the Pauli matrices 

= ( 1 ) ' = ( i 7 ) ' ^2 = ( -1 ) ' 

already have the right dimension and anticommutation properties, they can 
directly used to represent the Euclidean 7-matrices in D — 2 dimensions. 
One possible choice is 

70 = , 

71 = (7i . 

and for hermitian 75 

75 = ^7o7i = ■ 

In this representation the chiral projectors are diagonal. They are defined 

^«,L=^(1±75), (A.l) 

with the properties 

Pr,l-Pr,l, PlPr = PrPl = 0, (A.2) 

^,Pr,l = ±Pr,l, (A.3) 
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Pr + Pl = 1.^ Pr -~Pl = 75, (A.4) 
We also need the projectors defined using 70 

P± = l(l±7o). (A.5) 



These projectors become especially simple in another representation 

70 = 0-3 , 

71 = CTi . 

Then 

A. 1.2 Generators of SU(N) 

The Lie algebra su(A^) of SU{N) can be identified with the space of complex 
N X N matrices M satisfying 

M^ = M, Tr{M} = 0. (A.7) 

The generators of SU (N) can be identified with a basis A", a = 1, 2, ... , N^ — 
1 in this space normalised to 

Tr {X^X'j^CSab, C = 2. (A.8) 

They obey commutation relation (repeated indices are summed over) 

[A", A''] = 2i/"'''=A% (A.9) 

with the totally antisymmetric structure constants f"'^'^. 

Now every complex N x N matrix X can be expanded in the complete 
basis 

A^=|^1,A'^|, A^0,1,...,N'-1. (A.IO) 
The expansion reads 

X = I^A^ Tr{A^X} . (A.ll) 

A 



A.l Definitions 
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Expanding in such a way the matrix xf^'"' = 26ikSji one finds the identity 

(A")jj(A")ifc = 2Sik5ji - j^SijSik . (A. 12) 

Thus the quadratic Casimir operator is 

(A"A'^),, = C25.,, ^2= ^ ^ . (A.13) 

The following identities are usefull when evaluating flavour traces of Feyn- 
man diagrams 



Tr{A"A''} 


= NC2, 


(A.14) 


TrjA^A'^A^A''} 


= N{C,)\ 


(A.15) 


(A'^A^A'^)i, 


= (C2-2iV)(A%, 


(A.16) 


Tr{A"A^A"A''} 


= NC2{C2-2N), 


(A.17) 


Tr{A"A^} Tr{A"A^} 


= CNC2. 


(A.18) 



A. 1.3 Lattice notation 

Here we define the lattice difference operators. Via the factor A^ we are able 
to define the theory with general boundary conditions. (For details see j7j|.) 

dM^) = i[A;.V^(x + a/i) -^(x)] (A.19) 

= imx) - X-'^Pix - ai,)] (A.20) 

Xf, = e'^^^/^ , 00 = 0, -TT < 01 < TT , (A.21) 

p± = p,±9,/L (A.22) 

= isin(ap^) (A.23) 

= fsin(ap^/2), (A.24) 

The left action of difference operators is defined as 

^ix% = i[^(x + a/i)A;i-V^(x)] (A.25) 

H^)^; = i[?(x)-^(x-aA)A^]. (A.26) 
They are related to the right difference operators 

= -J2^{x)d;^Pix) (A.27) 

X X 

= -J2^ix)d,^ix). (A.28) 
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and can be used to define the left action of the Dirac operator 

D = W;+Kh,-oB;n,] (A.29) 

In particular, the action in the interior of the lattice can be written 

T—a L—a T—a L—a 

^ ^ '^{x){D + mo}^(x) ^o?^^ + mo}^(a;) . (A.31) 

xti=ax\=Q xo=axi=Q 

A. 2 Free theory 
A. 2.1 Formulae 

The positive energy plane wave solutions of the Dirac equation are 

V'(x) = e^P^ , Im po > , (A.32) 
with spatial momentum pi integer multiple of 27r/L in the range 

— 7r/a < pi < 7r/a , (A. 33) 

The energy 

Po—pQ— ^^{Pi) mod 27r/a , (A. 34) 

is constrained by Dirac equation giving 

(p+)2 + M(p+)2 = , (A.35) 

M{p) = mo + fp' , (A.36) 

which defines a;(gi) 

"'"'b"'«''1^2| l + a(mo + f@ / • 

This implies <^(Pi) > for mo > 0. 

The following amplitudes appear in the free propagator 

A{qi) = 1 + a(mo + ^qf) , (A.38) 
R{q) = M{q) {l - q-^^^'^'^^} - iqo {l + q-^^^'^^^^] . (A.39) 



A. 3 Four fermion operators 



105 



A. 3 Four fermion operators 
A. 3.1 Fierz transformation 

Fierz transformations connect products of Dirac bilinears by rearranging the 
order of the Dirac spinors. Thus some of these products are not independent. 
In two dimensions the set 

= {1,70,71,75} 

is normalised to 

tr(r^r^) = 25ab (A.40) 
and forms a complete basis for complex 2x2 matrices M^/j: 

^ A 

Considering four Dirac spinors '^j^, ■1/^2, and ip^ the general form of the 
Fierz identity is then 

(V^ir^V2)(V^3r^V'4) = c^cD(V^ir^V'4)(V^3r^V'2) ■ (a.4i) 

To fix the unknown coefficients Cq^ consider the left-hand side with all Dirac 
indices explicit 

= -iEra/3rf.r^«r^Atr(r^M)tr(r^ivr) 

C,D 

= -tE tr(r''r^r^r^)(^ir^^4)(V^3r^V^2) 

C,D 

where M^a = 'ipia'^Ma and M'g\ = ip^xi'W are 2x2 matrices, that are ex- 
panded in the third line in the basis given above. In the last line we used 

tr(r^M) = rJ,^i,V'4. = V^ir^^4 

and similar for tr(r^M'). The overall sign is due to the fact that the fields 
anticommute. Thus the coefficients are determined through 

^ -1 tr(r^r^r^r^). (a.42) 
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re 




1 


J-, 705 7i) 75 


1 

\ 


7/. 


1 

75 
lu 


2 

1 


75 


1,75 
lu 


~2 
1 

2 



Table A.l: Some coefRcients of the Fierz transformation Eq. l|A.4ip . 

We are mainly interested in the coefficients for = F^. One finds 

C^^oc(5c7D for all F^ , F^ . 

The 16 nonzero coefficients Cqq are collected in Table [AH 

For the subtraction of F^ = F^ = 1 and F^ = F'^ = 75 only the terms 
with 7^ survive 

(V^l^l)(^2^2) - (V^l75V^l)(V^2 75^2) = -X^(V^l7/.V^2) (V'2 7mV^i) • 

In the case of F"^ = 7^ the terms with F*" = 7^, cancel due to the implicit 
sum over /i and the fact that we have just two 7-matrices: 

(V^i7m^i)(V'27/.V'2) = -(V'i^2)(V'2'0i) + (V^i75V'2)(V'2 75''/'i) • 

A. 3. 2 Flavour mixing 

Using the expansion introduced in Eq. flA.lOl) the matrix Mjj = ip^Tipj, 
where F is a matrix contracting the Dirac indices of ip and ip, can be written 
as 

A 

and 

M,, M,, = - 5^ Tr { A^} F ^) F A^ ^) . 

i,j A,B 

Since the matrices A^ are normalised to 

Tr{A^A^} =25^5 
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one ends up with 



r v^,- V^^. r v^, = 1(V^ r v^f + i J] (V; r A« 

i,j a 
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Correlation functions 



B.l Properties of the free Wilson propagator 



B.1.1 Analytically known quantities 

We work in a half Fourier transformed space, e.g. the free propagator in the 
interior of the SF is 



a e*(«i--i) rnxMy)] = S{xo, yo, h) . (B.l) 



To evaluate the zeroth and first order diagrams we need the free propa- 
gators from boundary to boundary 



aJ2 e^^^(^^-^^) [C'(^i)C(yi)] = Hpi) P+ = K{p,) , (B.2) 

XI 

^iMm-x.) [((x.X'ivi)] = l,K\p,)j, , (B.3) 

XI 
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and from the boundaries to the bulk (and vice versa) 

«E e'^'^y^-^'^ mx)ayi)] = h{xo,pi) P+ (B.4) 

yi 

= H{xo,p,), (B.5) 

aJ2 e^"^(^^-^^) [CiyMx)] = l,H\xo,Pih, , (B.6) 
yi 

aJ2 e'"'^'-^'^ mx)C\yi)] = h{T-xo,Pi)P- (B.7) 

yi 

= H'{xo,Pi), (B.8) 

aY, e'"'^"'-^'^ [C{yi)i^{x)] = ^,H'\xo,PiH. (B.9) 
yi 

From the definition of the Wilson propagator (I4.4ip one infers 

k{h) = -2^k-^^e-^''t)T^ (B.IO) 

h{yo, ki) = hi{yo, ki) + 271/12(^/0, ki) , (B.ll) 
h{yo, k,) = M{k+) sinh(cu(A;+)(r - yo)) - tk^ cosh(cu(fc+)(T - yo)) , (B.12) 

h2{yo, ki) = -kt sinh(^(A;+)(T - yo)) • (B.13) 

For fci = and at zero bare mass we find 

fc(0) I , = , }r,rr.,r. + O(a^) , (B.14) 

V )\mo=o cosh{eT/L) ^ ^' ^ ' 

hiiyo, 0)L^^o = ^/^cosh(^(T - yo)/L) + 0{a') , (B.15) 

/^2(2/o,0)L^^o = 9/Lsmh{9{T - y,)/L) + 0{a^) , (B.16) 



2e-2a;(o+)r 



e/Lcosh{eT/L) 

mo=0 



^ +0(a'). (B.17) 



B.1.2 The propagator at zero distance 

Consider the free propagator 

S{x, y) = [^(x)^(y)] = {^P{x)lp{y)), , (B.18) 

which is one of the basic contractions in our computation. For coinciding 
arguments y ^ x and in the continuum theory, just by Euclidean invariance, 
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one infers that S{x, x) must be diagonal in Dirac space. This can be made 
more explicit by looking at the free propagator on the lattice in a periodic 
box. Since S{x, x) it is a local quantity, it is not sensitive to the special kind 
of boundary conditions, once the continuum limit has been taken. For an 
infinite periodic box the free propagator for Wilson fermions reads 

J-n/a 27r p^ + M{py 

where M{p) is defined as usual (cf. my notes). Except for — «7^p^, which 
is odd, and e''^^^~'^\ which is one for y — > x, all terms in this expression are 
even. Thus Spb{x,x) has no Dirac structure, i.e. is a multiple of the unit 
matrix in Dirac space, even at finite lattice spacing. This will change if we 
allow for general boundary conditions including a phase e*^. But it will still 
converge to the identity in the continuum limit. 

For the Schrodinger functional free propagator the same behaviour can 
be shown to hold. In particular this means for the dimensionless product 

a S{x, x) = ^Yl "^(^0' ^0' ^i) = ^(^o) + 0(«V^^ ^i) , (B.20) 

ki 

where B{xo) = Bi + B2{xo,6i) a/L + 0(a^/L^) is a c-number valued function 
and S{xo, i/q, ki) is the propagator in the half Fourier transformed space we 
do the actual calculation in. Its deviation from a diagonal form in Dirac 
space is of order a? / L"^ 6i. For the diagonal part we find 

B{xo) = 0.3849001 + [0.5000 + 0((T - 2xo)/L) + 0{9i)] a/L + 0(aV^^) , 

fB.21) 



B.1.3 Bubble reduction 

In the perturbative expansion the free propagator at zero distance appears 
sandwiched between two Dirac structures when two spinors of one and the 
same four fermion interaction are contracted and produce a bubble in the 
diagram (see also Fig. IB. II ) 

Tja [ip{x)^{x)~\Tj = a S{x, x) + Tj [aS{x,x),Ti] (B.22) 
= B{xo) + Fi{x) ayLHi + ... , (B.23) 

where on the left hand side the brackets denote a contraction while on the 
right hand side they denote a commutator. This means that in these cases 
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Q 



+ 0(aVL"9i) 



V 



Figure B.l: A four fermion interaction contracted to a bubble can be reduced to a 
insertion of the scalar density. 



the insertion of the four fermion interactions can be reduced to the insertion 
of the scalar density 

X 

a ^ (. . . {B{xo) + Fj{x) a^L' 6, + ...) ^{x) . . . >^ . (B.24) 



For ^1 = (periodic boundary conditions in space) this reduction is exact. 

The insertion of the scalar density in a correlation function in turn can 
be written as the derivation of this correlation function with respect to the 
mass parameter mo 



d 

darriQ 



0), = -aY,{Oi,{x)^{x))^ 



(B.25) 



B.2 Boundary-boundary correlation functions 
B.2.1 Free theory 

In a diagrammatic expansion we find at zeroth order one diagram for the 
two- and the four-point function (17.21) and (17.31) . In terms of the propagators 
introduced in the last section they read 

/f = tr,{K(0)} = fc(0), (B.26) 
/f = tr4ir(0)i^t(o)} =A;2(0). (B.27) 

Note that ff^ = (^fi^^^ ■ At zero bare mass we find 

+ 0(0^) . (B.28) 



AO) 

J2 



-0 cosh(^T/L) 
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B.2. 2 First order 

At first order there are two diagrams for and one diagram for /2 (cf. Fig. 
17.11) . Only diagram a is needed in the computation of the ratio (17. 5p . That 
is because with the help of Lemma [U diagram b can be written as f2^^ ■ f^^ 
and is hence canceled in the expansion. In terms of the above introduced 
propagators diagram a reads {H{yo) = H{yQ,0)) 

/i,a=E^^41' I = S,P,V, (B.29) 



filL=l E tTa{TBl,H\xoh5TiH'{xo)T'si,H'^ixoHTjH{xo)} . 

(B.30) 

The convention for coupling constants used in this appendix translates into 
the one of (I5.27p via 

cs = , cp = -5p, and cy = Qv ■ (B.31) 

Since {71, /i(?/o,Pi)} = in case of diagram a, it is easier to do the com- 
putation with V'p = V B = li (of course the result must be the same with 
r'^ = Ffi = 75, this can be used as a check). Then 



'2e-2c^{o+)T^ 



tr, {h{x^)P+^iP^h{xo)Tih{T - xo)P-iiP+h{T - xo)Ti} , (B.32) 

and 

tr, {. . . } = tr, {7ir,7i h{x^)P_h{x^) h{T - xo)P-h{T - Xq)} . (B.33) 
Using 

hP^h = hjP^ - hlP+ + i7i/ii/i2 , (B.34) 

this results in 

Tcj = 1 ^ iict{...}={hi{xQ)hi{T-XQ)-h2{xQ)h2{T-XQ)f 
Yp = tirf {...} = - {hi{xo)h2{T - xo) - /i2(xo)/ii(T - xojf 

Ty = 7o + 7i ^ tid {...} = - I {hi{xo)hi{T - xq) + h2{xo)h2{T - Xo) f 

+ {hi{xo)h2{T - Xq) + h2{xo)hi{T - xo)f I . 
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At zero bare mass we find 

rs = l ^ trd{...} =(^/L)^cosh2(^(T-2a;o)/L) + 0(a2) 
Tp = ^75^ tr, {...} = - smh.\e{T - 2xo)/L) + 0{a'') 

Ty = 70 + 71 ^ tr, {...} = - (6/1)' cosh{2eT/L) + 0{a^) . 

Taking the naive continuum limit of the sum in (IB.30I) for I = S,P we find 

T-a T 

J2 cosh^{9{T - 2y,)L) = - dyoCosh\9{T - 2y^)/L) + 0(a) 



a 
L 



yo=a 



1 sinh(2^T/L) + ^ + 0(g) , 



T-a „T 

I J2 smh\eiT - 2yo)L) = - dyo sinh2(^^(T - 2yo)/L) + 0(a) 
yo=a. ° 



1 sinh(2eT/L) - ^ + 0(a) 



Finally using flB.lTD 



(1) 

4,a 



2 - cp) A{d) + CS + CP- 2cvB{e)} + 0(a) 



(B.35) 

A{e) = ^ smH2eT/L) 1 , (B.36) 

B{e) = cosh(2^T/L) '^^^ 1 , (B.37) 
C{e) = cosh.\eT/L) 1 . (B.38) 



B.2.3 Second order 

At second order there are nine diagrams for /4 and three diagram for /2 
(cf . Fig. I7.2p . But because of Lemma [T] we need to evaluate only the five 
non-reducible diagrams a-e. 

The second order diagrams involve one momentum loop and are therefore 
treated numerically. For the logarithmic divergent part we find 



f(2) 
i4,a 



T ln(a/L) 



{{-4-cl + 2csCp)Bi9) 



- 4 4 B{e) + 4 cv{cs - cp) A{e) } + 0(1) , (B.39) 
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f(2) 
J4:,b 



T ln(a/L) 



mo=0 



|2 4 {A{e) + l)-2 4 {A{e) - l)-4cscp 



2L {Cie)f 271 

+ 4cv {cs {A{e) + 1) - CP {A{9) - 1)) } + 0(1) , (B.40) 



f(2) 
/4,c 



T ln(a/L) 



mo=0 



|(4 + 4 + 2c5Cp)S(^) 



2L ic{e)f 2% 

+ 44 -Acvics + cp)} + 0(1) , (B.41) 



f(2) 
^4,d 



-iVT ln(a/L) 



2L iC{e)f 2tx 



nio=0 



24 (A(^^) + l) 

-24 (A(^) - 1) } + 0(1), (B.42) 



f(0) 



= 0(1) (B.43) 

mo=0 



B.2. 4 Proof of Lemma [E 



A diagram /£j at order > in the expansion of /2 (cf. Eq. (17.21) ) takes the 
general form (the subscript i is only to label the diagram) 

/2? E ^''i{VnP+}, t = l,...,n, (B.44) 
where for n = 

Vo = KiO) , (B.45) 

and for n > 1 

Vn = h{T - yo) V{xl ■ k{ ■ It) hizo) , t=l,...,n, (B.46) 

2p-2<^(0+)T 

With /i(|/o) = ^(z/o, 0) as defined in Section IB.l.li For example, if n = 1 
there is only one possibility 

V, = h{T - xl) Ti, J2 ^(4' 4' ^1) T/i H4) ■ (B.48) 

ki 
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For arbitrary n> 1 the kernel V{xq ; k{ ; It) is a product of free propagators 
and 2n Dirac structures F/^ = {1, 75, 70 + 71} coming from the vertices. The 
free propagator is 75 hermitian 75S'(a;o, yo, ^1)^5 = S{yo,xo,ki). And since 
besides the propagators there is an even number of 7-matices in V, it also is 
75 hermitian. Then Vn has the property 

l5Vh5 = h{zo) V{4 ■k\-h)h{T-yo), t = n, . . . , 1 . (B.49) 

Using properties of the trace and that /2 is real we have the identity 

tU {Vn P+} = tU {75^175 P-] . (B.50) 

A reducible diagram at order n > takes the general form 



n . 



/i?=(f)" E t^d{rBP-DsP-r'^P+UrP+} ,t = l,...,n, r + i 

Xq ] ; It 

(B.51) 

Because of the projectors P± = |(1 ±70) only for boundary Dirac structures 
F^ , F^ G {71 , 75} Eq. (IB.51I1 does not vanish and reduces to 

/i? = (f)" E tT,{TsDXBP^UrP+} . (B.52) 

First we concentrate on proving the case n = 0. For r, s = we have 

Uo = K{0), 1^0 = 75^^^0)75. (B.53) 
and there is only one diagram 

/f = tT,{TBl5K\0)j,T'^P+K{0)P+} . (B.54) 
For Tb = ^'b we find 

/f = tr,{A^t(o)p^;^(o)P^}= tr4irt(0)P+}- tr,{ir(0)P+} = /f -/f , 

(B.55) 

where we used K^{0) = K{0) and which proves the Lemma in this case. For 
mixed choices for F^, F'^ there are factors ±i left which can be absorbed 
into Tb or F^. 

Now we turn to arbitrary n > 1. For r, s > 1 we find 

Ur = h{T - yo) U{xl ; kl ■ Q h{zo) , t; = 1, . . . , r (B.56) 
Vo^zo e {xl,...xl} (B.57) 
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and 

Ds = F^h{uo)D{x^;k';';I^)h{T-Vo), w = r + 1, . . . ,r + s (B.58) 
Wo,t^oe K+\...xS+^} (B.59) 

Because of the projectors P± in (1B.52I) and the fact that the Xq and k\- 
summations do not mix Ur and Dl the overall trace factorises into two inde- 
pendent traces 

. h.W .J . UV . J 

O' l' ^ O' l' '^ 

(B.60) 

The second factor is evidently an /2 diagram. Because of F'^ , e {71 , 75} 
and the cyclic properties of the trace the trace of the first factor can always 
be written as (one may have to absorb some factor ±i into F'^ or F^) 

tuiVBDsV'M = tu{DsP-} . (B.61) 

Comparing Eqs. (IR6T1) . ([RBSl) with flRBOl) . flR49l) also the first factor in 
Eq. (lOnli can be identified with an /2 diagram and we arrive at 

/i? = /g ■ & ■ (B.62) 

It should be obvious that this holds also for s > 1 , r = 0, or vice versa, which 
completes the proof of Lemma [H 



B.3 Boundary-to-interior correlation functions 
B.3.1 Free theory 

Here we consider tree level correlation functions of the type 

C2f,i(ri,FR;pi,xo) = -a^^ e^fi(^i-^i) (V^(x)Fi^(x) C(i/i)FrC(^i)>o , 

(B.63) 

where Fi is a matrix with Dirac and flavour indices specifying the quantum 
numbers of the current or density that is inserted at x. The disconnected 
part of Eq. (IB.63I 1 is proportional to 

iUj{S{x,x)Vi} tuj {P-SP+V^ - aP.^^if>tV^] (B.64) 

and the connected part 

iUj{P-STiSP+T^} . (B.65) 
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7i,75 
7iA",75A'^ 


7i,75 
7lA^75A'^ 



Table B.l: Nonvanishing Dirac and flavour structure for correlation functions with cur- 
rent/density insertions. 

Both expression vanish for Fr oc 1,7o. The disconnected part vanishes for 
nontrivial flavour structure of either Fr or Tj. The connected part can be 
written in terms of propagators between boundary and interior using (IB .41 - 

MM 

tidj {H\xo,pi)'y5TiH{xo,Pi)Tji'y5] . (B.66) 

The nonvanishing combinations of Fr, Fj are listed in the table Table IB.li 
Fi oc l,7o are excluded by parity. 
If we introduce the shorthand 

=aJ2 e'^^^y^-^^^ P^S{a, z,; a, y,)P+ , (B.67) 

2/1 

in the disconnected part the the correlation function Eq. (IB. 631 1 reads 

C2f,i(Fi, Fr; Pi, Xq) = CgJ - C^l , (B.68) 

CgU tTdj{H^{xo,Pi)^5TiHixo,Pi)Tj,^5} , (B.69) 
CSI = Yl ^^'i'f {^(^0' fci)ri} tr,,^ {L(pi)Fr} . (B.70) 

The half Fourier transformed of the propagator S is defined in Eq. (|5.56p . The 
minus sign infront of j is due to the additional trace of the disconnected 
diagram (cf. Section [5.5.21) . 
Since 

H{xo,Pi)'y5 = H{xo,pi)i'yoli = iH{xo,pi)'yi , (B-71) 

there are only two independent zero momentum correlation functions with 
current/density insertion and vanishing disconnected diagrams. Namely the 
insertion of the space component of the vector current 71 and pseudo-scalar 
density 75. Note that in two dimensions 7^75 = ie^,y7i/, therefore = iefj^iyV^ 
and fA = ifv, where /a is the correlation function of the insertion of the time 
component of the axial current. Comparing (IB. 631) with (16. 7p we have 

fA\^o) = ^<^2f,i(7o75,75;Pi,a;o) , (B.72) 
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and 



fp\xo) = ^C72f,i(75,75;pi,a;o) . 



(B.73) 



Using eqs. (14.421) . (IB.4p . (IB.67P and performing some algebra these correla- 
tion functions explicitly read 



(B.74) 



with 



(0,1) 



R{p+y 



^(Pl) A -2lo(p+)t\ "^l" 



^ 2ik+A{kl)R{k^ 



and 



(0,1) 



fi;^^^ = Npt§^{l-e-^-^^tW^j: 



R{p+) V ^ ) ^ 2A{kt)R{k^ 

|g-2a.(fc+)xo _ g-2i^(fe+){T-a;o)| _ 

Here we used the abbreviation 

M±(p+) = M{p+)±ip^ . 

Note that 



'id; + do)fr'i^o) = 2Mip^)cosh[aiuipt)] f^'^'ixo) 



(0,1), 



(B.75) 
(B.76) 

(B.77) 

(B.78) 
(B.79) 
(B.80) 

(B.81) 
(B.82) 



and 



E 



M{k+) cosh [acj(A;/ 



g-2Lj(fc+)a:o _ g-2<^{fc+){T-xo) 



(B.83) 
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Perturbation theory vs. Monte 
Carlo simulation 

C.l Full theory with bosonic auxiliary fields 

Consider the chiral Gross-Neveu model in two dimensions with Wilson fermi- 
ons and the interaction terms Oss ,Opp , Oyv- In the Schrodinger functional 
set up correlation functions are calculated from the generating functional Z 

(0) = ||oz|_ ^, (C.l) 

where 

Z = j D[^]D[^]exp{-5o-5, + (^,r/) + (rJ,^)} , (C.2) 

and 

In Eq. (IC.3P the interaction terms are defined only for < a;o < T. We 
could extent the definition to the boundary, but the fields there are no dy- 
namical variables and do not contribute to the following considerations. In 
fact the fermionic integration in Eq. (1C.2I) is only over fields ^(x), ipi^x) with 
< Xo < T. The free action 5*0 is defined in Eq. (14.341) with the Wilson- 
Dirac operator in presence of the boundaries (I4.35p . The choice of coupling 
constants translates into the one of (15.271) via 



gl = g'^ and gl = g^ - Sj, 



(C.4) 



122 



Perturbation theory vs. Monte Carlo simulation 



On the right hand side of Eq. ( IC.ll l the operators in O on the left hand side 
have been replaced by the corresponding functional derivatives. In pertur- 
bation theory one expands the interaction part of the generating functional 
in the coupling constant gj. In this way all sorts of vertices are added to the 
external operators in O. Here we want to introduce bosonic auxiliary fields 
absorbing the interaction terms and leaving a free fermionic theory with a 
modified Dirac operator. 

Introducing auxiliary fields defined by 

a = -glH' (C.5) 
= -gli'i'y^i/j (C.6) 
= -gl^j^ij (C.7) 

into the generating functional via 

1 oc / D[X]e-4(^-^^^^^-^)^ (C.8) 



the fermionic integration becomes Gaussian. 

But the generating functional now contains also the integration over the 
auxiliary fields and the free action now depends on them 

L—a 

So[a,n,Bf,] = ip{0,Xi)P^{a-fidii!{0,Xi) -i){a,xi)} 

xi=0 
T-a L-a 

+ ^ ^ i'{x){D + mo + a + ij^n + 7^5^}V^(x) 

L—a 

+ aJ2 ^{T, x,)P+{aj,dMT, Xi) - tP{T - a, Xi)}{.C.9) 

xi=0 

The the generating functional is up to constant factor 

Zoc j D[a]D[7i]D[B^,]Zo[a,7i,B^]e-'\ (C.IO) 
with the fermionic part of the generating functional 

Zo[a,7r,S^] = y"D[^]D[^] exp{-5o[(T,7r,5^] + (^,r7) + (7j,^)} , (C.ll) 
and the kinetic terms for the auxiliary fields 

xo=a a;i=0 



C.l Full theory with bosonic auxiliary fields 
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The operator O contains functional derivatives only with respect to fer- 
mionic fields. We can therefore write the expectation value Eq. (IC.lll 

(O) = I y D[a]D[7r]D[5^] (O), Zo[a, vr, 5^]^„,,=o e"^" , (C.13) 

where is taken with Zq^j^ti^B^ and the subscript a indicates, that it is 
still dependent on the auxiliary fields 

(0), = {OZo[a,7r,5^]}_^^„ . (C.14) 

Introducing the Dirac operator Da depending on the auxiliary fields 

= D + CT + 2757r + 7^5^ , (C.15) 

the fermionic functional integral ZQ[a^ vr, for vanishing source and bound- 
ary fields is the determinant 

Zo[(T, TT, 5^]p...^=o = (det(L'a + mo))^ . (C.16) 

The expectation value Eq. (IC.14P is just a sum of tree diagrams with the 
propagator Sa{x,y) defined through 

{Da + mo) Sa{x,y) = a~'^6^y , 0<Xo<T. (C.17) 

Using the contractions derived in Section every combination of field oper- 
ators in O can be expressed in terms of the propagator Sa{x, y). The remain- 
ing functional integral over the c-number valued auxiliary fields is accessible 
by numerical methods 

(O) = ^ y D[a]D[7r]D[fiJ (0)„ (det(D, + mo)f e"^" , (C.18) 

given the determinant can be well defined. Since we can choose a repre- 
sentation of the 7-matrices where 75 is imaginary and 70 , 71 real, Qxy = 
{Da + mo){x,y) is a real matrix. Therefore the eigenvalues have to come in 
complex conjugate pairs and the determinant is real 

detQ = (detQ)* = detQ^ (C.19) 

But because the determinant in Eq. (IC.18I) can be written as e'"^^'^'^*'^)'^, it 
has to be strictly positive. This, however, is the case for even N. 
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Perturbation theory vs. Monte Carlo simulation 



C.2 Results 

We calculate at finite lattice spacing 

am{xo) = 



2/p(xo) 



(C.20) 



in Monte Carlo simulations and in perturbation theory (PT). In the simu- 
lations we use the representation (1C.18I ) and a standard fermion algorithm 



42l |. Choosing 



Op = Qs 1 aAm = arriQ — arric 



with 



am, = -0.7698004(1) {N gl - g^) + 0{g^) 
we find in leading order perturbation theory (cf. Chapter [6l) 



am(xo) 



/in 



2fP 



2 ( ^1,J 



(1) 



^0 fpj 



+ 0{g' 



2/. 



(0) 



(C.21) 
(C.22) 

(C.23) 
(C.24) 



amo=aAm 



^2 ^0 damo fp 



(0)' 



,2/r 



2(/f )2 



(1) 



+ 



hi J hofpj 



2(/. 



(0)^2 



+ 0(/) 



(C.25) 



where 



ha = adoff\xo) , hi = adof^'jixo) , h2 = damoho , 



(1), 



(C.26) 



and in the second line we expanded arrio around aAm. 

For pi = (external momentum) and 9i = (for these values the discon- 
nected diagram vanishes) the tree level amplitudes can be computed analyt- 
ically, giving 

am{xo) = aM{p+) cosh(aw(p+)) + 0(/), . (C.27) 



Monte Carlo data and PT are plotted in Fig. IC.ll First order PT seems 
to be valid at g^ ^ 0.15. See [3j| for the same plot for a 12 x 13 lattice. 



C.2 Results 
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Cp = Cg, Cy = 0.01, aA m = 0.2 



LxT=6x7 

am. = A m + am* ' g 

MC: 10k updates 
PT: 0.1864- 0.1939 



0.05 0.1 0.15 0.2 0.25 0.3 

„2 



0.35 



Figure C.l: MC vs. PT. 
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